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FOREWORD 


I 

This  report  sunmarlzcs  the  results  of  the  vlndtunnel  test  of  an  os-  ] 

dilating  straked  wing  conducted  under  a  cooperative  program  of  research  ! 

I 

between  General  Dynamics  Fort  Worth  Division,  Fort  Worth,  Texas,  and  the  j 

National  Aerospace  Laboratory  (NLR) ,  The  Netherlands.  The  model  and  I 

1 

support  system  was  designed  and  fabricated  at  NLR  under  a  separate  pro¬ 
gram  with  General  Dynamics  and  NLR  funding.  The  test  preparation,  wind- 
tunnel  test  and  reporting  were  performed  at  NLR  under  Air  Force  Contract 
F33615-85-C-3013,  for  the  Flight  Dynamics  Laboratory  of  the  Air  Force 
Wright  Aeronautical  Laboratories,  Wrlght-Patterson  Air  Force  Base,  Ohio, 
the  work  was  administered  by  Mr.  D.W.  Kinsey  of  the  Aeromechanics  Divi¬ 
sion  (AFWAL/FIM) .  Additional  technical  monitoring  support  was  provided 
by  Mr.  T.  Cord  of  the  Flight  Control  Division  (AFWAL/FIG) . 

The  program  manager  and  principal  investigator  was  Dr.  A.M.  Cunningham 
Jr.  at  General  Dynamics  and  Mr.  R.G.  den  Boer  was  the  principal  investi¬ 
gator  at  NLR.  Mr.  den  Boer  was  assisted  by  the  following  NLR  specialists: 

C.S.G.  Dogger,  E.G.M.  Geurts,  A.J.  Persoon,  A.P.  Retll  and  R.J.  Zwaan. 

This  report  consists  of  six  parts.  Part  I  presents  a  general  description 
of  the  model  and  test  program  and  a  discussion  of  the  results.  Part  II 
contains  the  steady  pressure  distribution  plots  and  the  major  part  of 
the  zeroth  and  first  order  harmonic  unsteady  pressure  distribution  plots. 

Part  III  contains  the  remainder  of  the  unsteady  pressure  distribution 
plots  and  plots  of  the  steady  and  the  zeroth  and  first  order  harmonic 
unsteady  overall  loads.  Part  IV  contains  time  history  plots  of  the  un¬ 
steady  pressures  and  overall  loads.  Part  V  contains  power  spectral  den¬ 
sity  plots  of  the  overall  loads  at  harmonic  oscillation  and  time  history 
plots  of  overall  loads  for  (1-cos)  model  motions.  Part  VI  contains 
results  of  the  flow  visualization  program. 
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Fig.  45  Influence  of  frequency  on  the  time  history  of  the 
vortex  core  positions 

Fig.  46  The  effect  of  the  wire  suspension  blocks 
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1  INTRODUCTION 

l.I  Motive 

Straked  vlngs  are  currently  beconlng  common  features  of  advanced 
fighter  aircraft.  The  strakes  are  designed  to  generate  vortices  from  high¬ 
ly  swept  leading  edges  which  stabilize  flow  over  the  wing  and  provide 
additional  lift  up  to  high  angles  of  attack.  In  this  way  the  strakes  con¬ 
tribute  much  to  high  maneuvrabllity.  The  vortex  lift  capability  of  straked 
wings  has  been  extensively  explored  and  experimental  data  concerning  aero- 
d3niamlc  loading  are  available  for  various  planforms  and  Mach  numbers. 
Advanced  calculation  methods  to  predict  the  vortex  flow  are  maturing. 

The  knowledge  of  unsteady  loading  on  straked  wings  Is  less  developed, 
both  In  the  cases  where  the  loading  Is  due  to  wing  oscillations  -  as 
required  for  aircraft  stability  and  flutter  analysis  -  and  In  cases  where 
fluctuations  In  the  flow  are  Induced  by  vortex  burst  (or  vortex  breakdown) 
as  required  for  stall  and  buffet  predictions.  The  common  practice  In  flut¬ 
ter  analysis,  for  example,  Is  that  flutter  clearance  calculations  for 
straked  wing  fighters  are  still  based  on  attached  flow  without  leading 
edge  vortices.  This  Is  necessary  because  a  lack  of  adequate  calculation 
methods  for  unsteady  vortex  flow  which  must  be  validated  by  experiment. 

The  Immediate  question  then  arises  If  such  a  procedure  can  guarantee  suf¬ 
ficiently  conservative  predictions,  or.  If  so,  Is  It  unreasonably  conser¬ 
vative  . 

This  situation  has  given  rise  to  a  cooperative  project  of  General 
Dynamics  (GD)  and  the  National  Aerospace  Laboratory  (NLR) ,  with  funding 
provided  by  the  Air  Force  Wright  Aeronautical  Laboratories  (AFWAL) ,  con¬ 
cerning  a  wlndtunnel  experiment  with  an  oscillating  straked  wing.  This 
experiment  will  extend  the  knowledge  of  unsteady  vortex  flow,  and  provide 
a  data  base  for  the  validation  of  computer  codes  for  unsteady  airloads. 


1.2  Some  physical  aspects 

A  brief  description  of  some  physical  aspects  of  the  unsteady  vortex 
flow  Is  given  here.  First  the  main  characteristics  of  the  steady  flow  are 
considered.  In  figure  2  a  straked  wing  at  Incidence  Is  shown.  Vortices 
are  shed  from  the  leading  edges  of  the  strake  and  the  wing.  The  sharp 
strake  leading  edges  generate  vortex  sheets  even  at  low  Incidence  which 
roll  up  spirally  Into  the  strake  vortices  and  flow  downstream  over  the 
wing.  The  vortices  Induce  strong  lateral  velocities  at  the  strake  upper 
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surface,  giving  rlaa  to  suction  peaks  at  the  position  of  the  vortex  cores 
(see  upper  left  Inset).  When  the  lateral  velocities  are  large  enough, 
secondary  flow  separations  occur,  leading  to  secondary  vortices 
spiralling  opposite  to  the  primary  vortices. 

At  moderate  Incidences  vortex  sheets  start  to  develop  from  the  wing 
leading  edges,  starting  at  the  kinks.  As  a  result,  the  spanvlse  pressure 
distribution  now  shows  two  pressure  peaks  (see  upper  right  inset) . 

At  higher  Incidences  vortex  burst  or  vortex  breakdown  occurs  Initial¬ 
ly  for  the  wing  vortices  followed  later  by  the  strake  vortices.  An  impor¬ 
tant  consequence  of  vortex  burst  Is  that  the  corresponding  suction  peaks 
become  weaker  and  that  the  vortices  lose  their  ability  to  produce  addi¬ 
tional  lift.  A  normal  behaviour  of  vortex  burst  Is  that  It  will  move  up¬ 
stream  when  the  Incidence  Increases.  At  still  higher  Incidences  large- 
scale  boundary  layer  or  stall  separation  occurs,  starting  often  at  the 
trailing  edge. 

The  explanation  of  the  above  vortex  flow  becomes  increasingly  com¬ 
plicated  In  case  of  Interactions  of  strake  and  wing  vortices,  their  in¬ 
fluence  on  vortex  burst  and  flow  separation  and,  at  high  enough  speeds 
their  interaction  with  shock  waves.  The  analysis  of  these  phenomena  is 
still  the  subject  of  many  Investigations  (e.g.  Ref.  1). 

When  the  straked  delta  wing  Is  oscillating,  the  strength  and  the 
position  of  the  wing  and  strake  vortices  will  oscillate.  As  the  vortices 
are  being  fed  through  the  vortex  sheets  emanating  from  the  leading  edges, 
it  is  to  be  expected  that  the  oscillations  of  vortex  strength  and  position 
will  lag  the  wing  oscillation.  Phenomena  like  vortex  burst  and  stall  sepa¬ 
ration  have  shown  hysteresis  effects  in  steady  measurements  so  that  It 
might  be  conjectured  that  in  the  unsteady  case  these  effects  will  cause 
an  additional  lagging. 


1.3  Objectives 

With  the  description  of  the  problem  area  In  section  1.1  and  the  vari¬ 
ous  physical  aspects  In  section  1.2  as  background  the  objectives  of  the 
windtunnel  experiment  were  formulated  as  follows: 

-  to  obtain  a  physical  Insight  of  the  vortex  flow 

-  to  set  up  a  data  base  of  unsteady  aerodynamic  data  for  computer  code 
validation;  and 

-  to  study  the  djmamlc  aspects  of  vortex  bursts,  up  to  high  angles  of 
attack  at  which  vortex  breakdown  occurs  close  to  the  strake  leading 
edge. 
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2.1 


2.2 


In  order  to  achieve  these  objectives!  sone  basic  decisions  were  made: 

a.  The  experiment  was  restricted  to  low  speeds  to  keep  the  already  compli¬ 
cated  flow  from  being  further  complicated  by  compressibility  effects. 
Continued  experiments  In  future  may  hopefully  cover  the  transonic  speed 
range . 

b.  Preference  was  given  to  a  balance  to  measure  the  overall  aerodynamic 
loads.  Because  of  the  large  pressure  gradients  It  was  expected  that  no 
practical  grid  of  pressure  pickups  Installed  In  the  wlndtunnel  model 
would  be  dense  enough  to  provide  accurate  aerodynamic  loads  through 
Integration.  Measurements  of  pressure  distributions  would  only  be 
applied  In  a  few  strategically  chosen  sections. 

c.  The  aerodynamic  load  and  pressure  measurements  would  be  supplemented 
with  Information  about  the  flow  above  the  model  upper  side.  Flow  visua¬ 
lization  was  used  to  study  the  position  of  the  vortices  and  vortex 
cores,  development  of  vortex  burst,  etc.  Application  of  a  laser  light 
screen  technique  In  a  previous  wlndtunnel  test  had  already  proven  to 

be  feasible  (see  Ref.  2)  to  study  the  position  of  the  vortex  cores  in 
steady  flow  (see  Fig.  3) .  Also  the  means  to  apply  a  chopped  laser 
light  screen  were  Investigated  for  unsteady  flow. 


TEST  SETUP 


Wlndtunnel 

The  tests  were  conducted  In  the  low  speed  wlndtunnel  L5T  3x2.25  m^ 
of  NLR  In  the  Moordoostpolder  (North-East  Polder)  of  The  Netherlands.  The 
tunnel  Is  a  closed-type  wlndtunnel  and  has  Interchangeable  test  sections 
with  a  total  length  of  8.75  m.  The  test  sections  are  provided  with  turn¬ 
tables  In  the  floor  and  celling.  The  tunnel  has  a  maximum  attainable  wind- 
speed  of  about  85  m/s. 


Model  and  model  support 

An  overall  view  of  the  test  setup  Is  shown  In  figure  4.  The  geometry 
of  the  model  Is  Indicated  In  figure  5.  It  Is  a  straked  delta  wing,  with 
a  span  of  800  mm,  a  root  chord  of  785.5  mm  and  tip  chords  of  120  mm.  As 
the  model  was  not  Intended  to  represent  a  real  aircraft,  the  parts  which 
might  complicate  the  aerodynamics,  like  fuselage  and  empennage,  were  left 
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out.  The  outboard  wing  has  an  NACA  64A005  airfoil  and  a  leading  edge  sweep 
of  40  deg.  The  strake  has  a  sharp  leading  edge  with  a  sweep  of  76  deg. 

The  spanwise  cross-section  of  the  strake  has  a  diamond  shape  with  a  half 
top  angle  of  11.4  deg  (see  figure  5).  At  the  centre,  the  diamond  shape  is 
rounded  with  a  radius  which  varies  linearly  from  0  mm  at  the  apex  x-0  mm, 
to  a  radius  of  80  an  at  xa333.7  mm  and  back  to  a  radius  of  0  mm  at 
x*437.0  mm.  At  the  kink  in  the  leading  edge,  where  the  strake  joins  the 
wing,  the  thickness  distribution  is  smoothed. 

The  middle  of  the  main  wing  was  thickened  to  accommodate  the  balance.  The 
sides  of  this  thicker  region  were  rounded  with  a  radius  of  80  mm  and  the 
top  is  a  flat  surface,  parallel  to  the  model  reference  plane  (see  Fig. 5). 

A  more  detailed  description  of  the  model  geometry  is  presented  in  appen¬ 
dix  A. 

The  model  was  designed  (see  Ref.  3)  and  fabricated  at  NLR.  It  was  made  of 
a  magnesium  alloy  to  reduce  inertia  loads.  The  overall  mass  of  the  model. 
Including  its  Instrumentation,  but  without  the  balance  was  5.9  kg.  The 
instrumentation  is  described  in  section  2.3.  For  the  visualization  tests 
a  smoke  tube  was  attached  to  the  lower  side  of  the  model,  with  its  opening 
very  close  to  the  apex  (see  Fig.  Sb-c) .  The  optimal  position  of  the  smoke 
tube  was  determined  in  a  separate  test  (see  Ref.  4),  although  during  the 
final  test  the  tube  could  be  shifted  in  x-  and  z-dlrectlons  and  its  inci¬ 
dence  could  be  changed. 

The  support  mechanism  (see  Figs.  4,6a),  also  designed  and  built  at 
NLR,  was  mounted  on  the  turntable  in  the  floor  of  the  test  section  to 
allow  model  sideslip.  By  a  system  of  struts  the  model  was  supported  in 
such  a  way  that  it  could  perform  a  pitching  motion  about  an  axis  at 
73.3  X  root  chord.  The  support  elements  were  aerodynamlcally  shaped  and 
two  of  them  were  used  as  conductors  for  the  instrumentation  wiring. 


3  Model  instrumentation 

The  instrumentation  of  the  model  consisted  of: 

-  1  six-component  balance, 

1  displacement  transducer, 

I  temperature  transducer, 

-  9  accelerometers, 

-  42  miniature  pressure  transducers. 

The  location  of  the  instrumentation  in  the  model  is  indicated  In  table  1 
and  figure  5. 


For  measuring  the  forces  and  moments  an  existing  six-component 
balance  (NLR  635)  was  used.  With  this  balance  forces  and  moments  can  be 
measured  up  to:  450  N  In  tangential  force,  1200  N  In  side  force,  3000  N 
In  normal  force,  90  Nm  In  roll,  110  Mm  In  pitch  and  70  Mm  In  yaw. 

With  a  Sangamo  AFG  5.0  S  linear  variable  differential  transducer 
(LVDT)  mounted  between  the  model  and  support,  the  oscillation  amplitude 
and  the  mean  steady  wing  Incidence  were  measured.  This  provided  the  cor¬ 
rect  geometric  Incidence  which  Included  the  deformation  of  the  balance. 

The  vibration  modes  were  measured  with  five  Endevco  2220  C  accelero¬ 
meters  of  the  piezo-electric  type  (range:  up  to  10E4  m/s^)  and  with  four 
Kullte  GY-155  accelerometers  of  the  plezo-reslstive  type  (range:  up  to 
500  m/s^) . 

Three  types  of  pressure  transducers  were  used: 

-  10  Endevco  8507-5 

-  12  Kullte  CQL-080-5D 

-  20  Kullte  XCS-093-5D 

These  transducers  were  mounted  In  such  a  way  that  they  were  electrically 
Isolated,  free  of  local  model  deformation  and  not  Influenced  by  the  model 
accelerations.  They  were  divided  over  four  strategically  chosen  sections: 

-  a  spanwlse  section  on  the  strake  to  obtain  data  In  the  conical  flow 
region  for  the  verification  of  conical  flow  phenomena  (10  pressure 
transducers  at  x/cr  *  0.4042), 

-  a  spanwlse  section  just  behind  the  kink  to  show  the  development  of  the 
leading  edge  wing  vortex  starting  from  the  kink  (18  pressure  transducers 
at  x/cr  -  0.6588) , 

-  a  spanwlse  section  at  the  rear  part  of  the  main  wing  for  measurement  of 
buffet  phenomena  (8  pressure  transducers  at  x/cr  -  0.9682), 

-  a  chordwlse  section  at  the  main  wing  panel  to  show  the  development  of 
vortex  burst  as  a  function  of  Incidence  (8  pressure  transducers  at 
2y/b  -  0.4000). 

The  sensitivity  of  the  pressure  transducers  shows  a  small  variation  with 
temperature.  By  measuring  the  temperature  of  the  model  with  a  Unlcurve 
thermistor,  UUT  45J1,  the  correct  sensitivity  of  the  pressure  transducers 
could  be  applied. 


Model  excitation 

Excitation  was  provided  by  an  electro-hydraulic  shaker  system, 
designed  by  Keelavlte  Hydraulics  Ltd.  It  consists  of  a  hydraulic  power 
supply,  a  combined  linear  actuator  and  servovalve  and  a  feedback  control 
unit  (see  Ref.  5).  The  hydraulic  actuator  could  deliver  a  maximum  static 
force  of  13000  N  and  a  dynamic  force  of  8000  N  at  a  total  piston  stroke 

of  35  mm  at  low  frequencies  to  16  mm  at  16  Hz. 

The  hydraulic  actuator  was  suspended  In  a  box  which  was  bolted  rigid¬ 
ly  to  the  turntable.  The  piston  was  connected  to  a  crank  which  converted 
the  linear  motion  Into  a  rotational  motion.  By  the  driving  rod  this  motion 
was  transmitted  to  the  yoke,  which  supported  one  side  of  the  balance.  On 
the  other  side  of  the  balance  the  model  was  clamped.  During  most  of  the 

runs  the  model  was  forced  Into  a  sinusoidal  motion.  In  a  limited  number 

of  runs,  a  (1-cos)  motion  was  used. 

Equipment  for  measurement  of  overall  loads  and  pressures 

The  wlndtunnel  tests  were  performed  using  a  computer  controlled 
multi-channel  transfer  function  analyzer,  called  PHAROS  (Processor  for 
Harmonic  Analysis  of  the  Response  of  Oscillating  Surfaces).  A  description 
of  this  system  Is  given  In  reference  6.  The  system  Is  capable  of  analyzing 
the  Incoming  data  from  48  channels  simultaneously.  By  means  of  a  switch 
panel  this  process  was  conducted  twice  at  each  test  point  so  that  96 
different  signals  could  be  recorded. 

In  figure  7  a  block-diagram  of  the  test  setup  is  shown.  Including 
the  PHAROS  system.  The  two-phase  oscillator  of  PHAROS  controlled  the 
hydraulic  actuator,  which  provided  the  model  excitation.  The  response  sig¬ 
nals  of  the  Instrumentation  were  acquired  by  PHAROS  for  analysis.  The 
electrical  signals  were  fed  through  conditioners  Into  transfer  function 
analyzers  to  obtain  the  mean  component  and  the  real  and  Imaginary  parts 
of  the  harmonic  components.  These  data  were  then  stored  on  a  disc  of  the 
computer  and  a  quick- look  analysis  was  made.  The  analysis  of  each  test  run 
was  performed  In  about  3  minutes.  Including  plots  and  tables  of  all 
measured  quantities  (see  also  section  4.1).  Thus,  Immediate  access  to 
detailed  preliminary  pressure  and  force  data  was  possible  during  the 


Equlpmtnt  for  flow  visualization 

In  figure  8  a  schematic  overview  is  given  of  the  visualization  test 
setup.  On  the  left-hand  side  a  top  view  is  shown  of  the  wlndtunnel  test 
section  with  the  wlndtunnel  model  and  the  flow  visualization  equipment 
beside  it.  The  smoke  tube  underneath  the  strake  injected  smoke  into  the 
flow  in  upstream  direction  near  the  stagnation  point.  In  this  way  the 
smoke,  mixed  with  air,  was  sucked  into  the  vortices  over  the  model.  By 
means  of  a  5  Watt  argon  ion  laser  and  a  cylindrical  lens,  a  light  screen 
was  formed  which  was  perpendicular  to  the  model  reference  plane  when  the 
model  was  placed  at  its  mean  incidence.  The  flow  patterns  were  made 
visible  by  the  light  scattered  by  the  smoke  particles  in  the  light  sheet 
which  is  shown  in  the  figure  at  the  right. 

In  order  to  record  the  flow  characteristics  at  fixed  phase  angles 
with  respect  to  the  model  motion,  a  chopped  laser  light  screen  was  used. 

By  means  of  an  acousto-optic  modulator  the  laser  light  screen  was  made 
Intermlttant .  The  modulator  was  controlled  by  a  special  device  developed 
at  NLR,  which  generated  electronic  pulses  in  phase  with  the  same  digital 
oscillator  signal  that  controlled  the  model  motion.  The  pulse  signal  was 
then  converted  into  an  amplitude-modulated  high  frequency  signal  by  a 
driver  which  excited  the  acousto-optic  modulator.  When  no  modulation  was 
applied,  the  laser  beam  was  Interrupted  by  a  mirror  and  reflected  to  a 
black  absorber.  In  case  of  lOOZ  modulation  tne  modulator  deflected  the 
laser  beam  about  6  mRad.  After  that  the  beam  passed  the  Interrupt  mirror 
and  was  led  into  the  optical  parts,  including  a  cylindrical  lens,  to  pro¬ 
duce  the  light  screen. 

The  optical  elements  could  be  rotated  about  the  optical  axis  to  place 
the  light  screen  perpendicular  to  the  wing,  when  the  model  was  at  mean 
incidence.  At  a  distance  of  2.50  m  downstream  of  the  rotation  axis  of  the 
model  an  Olympus  0H4  (3Smm  SLR)  photocamera  was  Installed  with  its  optical 
axis  in  the  symmetry  plane  at  the  same  level  as  the  rotation  axis  of  the 
model.  Photographs  were  taken  by  remote  control  with  a  250  fllmback  and  a 
100  mm  lens.  The  photographs  were  taken  to  determine  the  vortex  core  posi¬ 
tion  at  different  phase  angles  of  the  model  motion. 

To  enable  the  study  of  phenomena  like  vortex  burst,  also  video  record¬ 
ings  were  made  with  a  Charged  Coupled  Device  (CCD)  camera,  positioned 
beside  the  tunnel  wall. 
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PREPARATORY  TESTS 

The  model  was  supported  by  a  structure  consisting  of  several  struts 
(see  Flg>  4).  To  estimate  the  Interference  with  the  flow,  steady  measure¬ 
ments  were  performed  In  a  separate  preparatory  test  with  the  model  sus¬ 
pended  by  wires  to  an  overhead  balance  system  (see  Ref.  7).  Support  Inter¬ 
ference  corrections  were  derived,  which  were  applied  later  In  postproces¬ 
sing  the  results  of  the  main  test.  In  the  results  of  this  preparatory  test 
as  presented  In  reference  7,  the  pitching  moment  Is  affected  by  the  drag 
of  the  wires,  which  yield  an  additional  pitching  moment.  Only  the  force 
In  the  wire  attached  to  the  front  of  the  model  Is  used  to  derive  the  cor¬ 
rect  pitching  moment  as  presented  In  appendix  C. 

In  a  separate  test  (Ref.  8),  the  vibration  modes  of  the  model  (on 
the  balance)  and  the  support  system  were  measured.  All  resonance  frequen¬ 
cies  were  far  beyond  the  highest  excitation  frequency  in  the  teat  program 
(see  table  2).  Therefore,  the  influence  of  elastic  deformations  of  model 
and  support  mechanism  was  negligible. 

Another  preparatory  teat  was  carried  out  to  determine  the  optimal 
position  of  the  smoke  tube  with  respect  to  the  model.  Use  was  made  of  an 
oscillating  dummy  model  with  an  adjustable  smoke  tube.  A  steady  laser 
light  screen  was  applied.  The  test  was  performed  In  a  small  wlndtunnel  at 
the  Delft  University  of  Technology,  Department  of  Aerospace  Engineering. 

In  this  test  the  optimal  conditions  for  the  smoke  generation  were  also 
Investigated  (see  Ref.  4).  The  same  dummy  model  was  used  In  a  later  wind- 
tunnel  test  In  which  the  equipment  and  procedures  for  the  unsteady  flow 
visualization  were  tested. 


4  PROCEDURES 

4.1  Measurement  of  overall  loads  and  pressure  distributions 

4.1.1  Transfer  functions 

The  main  objective  of  these  measurements  was  to  determine  the  trans¬ 
fer  functions  between  the  mechanical  motion  of  the  model  as  Input  and  the 
pressures  and  overall  loads  as  output.  By  using  the  PHAROS  system  (see 
section  2.5),  the  zeroth  and  first  harmonics  of  the  measured  signals  were 
obtained  and  stored  on  the  discs  of  the  PHAROS  computer.  Next,  the  data 
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of  the  balance  were  corrected  for  Inertial  and  gravitational  loads, 
according  to  the  method  of  appendix  B.  All  quantities  were  normalized 
with  the  response  of  the  displacement  transducer  for  the  model  motion. 

Then  a  quick-look  printout  of  all  measured  quantities  was  produced  on¬ 
line  In  tables  and  plots.  The  complete  procedure  of  data  acquisition  and 
presentation  required  about  3  minutes  for  each  test  run. 

After  the  wlndtunnel  test  the  final  data  reduction  of  these  harmonic 
data  was  made  on  the  HP21/MXE  computer  of  the  PHAROS  system.  The  aerody¬ 
namic  quantities  obtained  after  each  test  run  were: 

-  one  chordwlse  and  three  spanwlse  distributions  of  the  mean  steady 
pressure  coefficients  (Cp)m 

-  one  chordwlse  and  three  spanwlse  distributions  of  the  unsteady  pressure 
coefficients  Re(Cp)  and  Im(Cp),  normalized  with  the  angular  displacement 
of  the  wing  derived  from  the  output  of  the  LVDT 

-  zeroth  and  first  harmonics  of  force  and  moment  coefficients,  measured 
with  the  NLR  635  balance,  again  normalized  as  mentioned  above 

-  amplitudes  and  displacements,  derived  from  accelerometer  signals  norma¬ 
lized  as  mentioned  above. 

4.1.2  Time  recordings 

Time  recordings  of  pressures  and  overall  loads  were  also  made  to 
enable  the  study  of: 

-  higher  harmonics  In  case  of  strong  nonllnearltles 

-  power  spectra  In  case  of  vortex  burst  and  strong  flow  separation 
effects 

-  cross-correlation  functions  In  space  and  time  when  following  the  develop¬ 
ment  of  vortex  burst  and  boundary  layer  separation. 

In  the  processing  of  the  recordings  of  the  overall  loads  to  time  history 
plots  of  overall  aerodynamic  coefficients,  the  procedure  Indicated  In 
appendix  B  Is  used. 


4.2  Incidence,  blockage  and  dynamic  pressure  distribution 


4.2.1  Incidence 


In  this  wlndtunnel  test,  a  very  large  range  of  Incidences  and  ampli¬ 
tudes  was  covered.  Although  the  output  of  the  LVDT  was  proportional  with 
Its  elongation,  the  elongation  Itself  varied  nonllnearly  with  the  Incl- 
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dence.  During  Che  experiment  this  nonlinear  relationship  was  not  yet 
determined  and  a  preliminary  relationship  between  output  of  the  LVDT  and 
incidence  and  amplitude  was  used.  For  chat  reason  these  preliminary  ad¬ 
justed  values  (aa  and  daa  are  presented  in  the  tables  which  give  an  over¬ 
view  of  the  test  program  (tables  4  to  12).  However,  in  the  postprocessing 
the  nonlinearity  was  taken  into  account  and  the  correct  geometric  inci¬ 
dence  (ag)  and  amplitude  (do)  were  obtained.  These  values  (measured  by 
the  LVDT)  were  not  Influenced  by  the  deformation  of  the  balance.  Next, 
additional  corrections  to  og  were  applied  to  enable  comparison  to  free- 
f light  conditions,  as  the  incidence  was  affected  by  the  presence  of  the 
support  system  and  the  windCunnel  walls.  First  a  zero-lift  correction 
was  applied  to  correct  for  the  Influence  of  the  support  system.  This  cor¬ 
rection  (-0.05  deg)  la  the  difference  in  incidence  at  zero  lift,  between 
the  steady  tests  in  wires  (see  Ref.  7)  and  the  steady  tests  on  the  strut 
support.  In  both  testa  the  wing  model  was  equipped  with  wire  suspension 
blocks.  Secondly  an  upwash  correction  according  the  method  of  references 
9  and  10  was  applied  to  take  into  account  the  influence  of  the  tunnel 
walls . 


4.2.2  Blockage 

Due  to  the  presence  of  the  model  in  the  test  section,  the  effective 
dynamic  pressure  is  increased  by  so-called  "blockage"  effects.  The  solid 
blockage  can  be  neglected,  due  to  the  small  volume  of  the  model.  The  wake 
blockage  is  primarily  caused  by  the  flow  separation  (bluff  body  drag)  and 
can  be  estimated  from  the  lift  vs.  drag  curve  according  to  the  method  of 
references  9,  10  and  11. 

During  the  windtunnel  test  the  lift  vs.  drag  curve,  determined  pre¬ 
viously  for  the  model  suspended  by  wires  (see  Ref.  7),  was  used  to  adjust 
the  tunnel  speed  in  such  a  way  that  the  dynamic  pressure,  corrected  for 
blockage  effects,  was  almost  independent  of  the  incidence.  In  the  postpro¬ 
cessing,  correction  for  the  blockage  effect  is  done  on  basis  of  the  lift 
vs.  drag  curve  as  measured  during  the  steady  part  of  the  test  with  the 
model  on  the  struts. 


A. 2. 3  Dynamic  pressure  distribution 

In  order  to  get  some  Information  of  the  dynamic  pressure  distribution 
around  the  model  In  relation  to  the  dynamic  pressure,  measured  by  the 
tunnel  reference  system,  some  dynamic  pressures  were  measured  for  zero- 
lift  conditions.  The  pressures  were  measured  In  a  plane  normal  to  the 
longitudinal  test  section  axis  (y'-z*  plane)  at  the  position  of  the  pitch¬ 
ing  axis  of  the  model.  The  pressures  measured  by  the  tunnel  reference 
system  were  corrected  to  values  in  the  empty  test  section.  The  djrnamlc 
pressures  measured  by  the  pitot  static  tube,  divided  by  the  dynamic  pres¬ 
sure  In  the  empty  teat  section  (Qref)  are  given  In  table  3  and  figure  9. 
These  data  were  Intended  for  possible  future  Investigations,  but  were  not 
used  In  this  report. 

4.3  Flow  visualization 

In  figure  10  the  principle  of  the  chopped  laser  light  screen  tech¬ 
nique  Is  shown.  The  upper  part  of  the  figure  shows  the  time  history  of 
the  model  motion  and  the  lower  part  shows  the  time  history  of  the  light 
pulses  which  generate  the  laser  light  screen.  The  phase  difference 
between  model  motion  and  light  pulses  could  be  varied  over  the  entire 
period  of  360  deg.  The  pulse  duration  could  be  varied  between  1  and 
8  deg. 

For  recording  vortex  position  data,  the  Olympus  OM  4  still  photo 
camera  was  used  (see  figure  8)  with  Kodak  Trl-X  as  film,  upgraded  to  1600 
ASA. 

As  the  light  screen  was  present  only  during  a  fraction  of  the  cycle 
(2.2  Z  maximum),  several  pulses  were  needed  for  one  exposure.  Consequently 
the  quality  of  the  photographs  was  strongly  Influenced  by  the  model  motion 
during  the  light  pulses  (depending  on  pitch  rate  and  pulse  width) ,  and  the 
time-averaging  effect,  l.e.  the  number  of  pulses  needed  for  one  exposure. 
This  was  not  a  serious  drawback.  In  a  preparatory  test  (see  section  3) , 
the  photographs  taken  were  very  suitable  for  the  purpose  of  determining 
vortex  core  positions  at  different  phase  angles  of  the  model  motion. 

For  qualitative  dynamic  Information  the  CCD  video  camera  was  used. 
Although  this  camera  has  a  smaller  resolution  than  the  photocamera,  a 
complete  picture  could  be  generated  from  each  light  pulse,  due  to  the 
high  sensitivity. 


, **  -  t 


-25- 

5  TEST  PROGRAM 

5.1  Overall  loads  and  pressure  measurements 

5.1.1  Steady  tests 

The  steady  tests,  performed  at  zero  sideslip  and  at  a  velocity  of 
80  m/s,  covered  an  angle  of  attack  range  of  -8  to  50  deg.  To  enable  com¬ 
parison  with  the  steady  test  of  the  model  suspended  by  wires  (see  Ref.  7), 
steady  tests  were  performed  with  the  model  supported  by  the  struts,  both 
with  and  without  wire  suspension  blocks  mounted  In  the  wings.  The  effect 
of  the  wire  suspension  blocks  had  to  be  examined  to  determine  their 
effect  on  Che  zero  lift  correction  (see  section  4.2.1). 

A  survey  of  the  steady  test  program,  Including  the  run  numbers,  Is 
shown  In  Cable  4. 


5.1.2  Unsteady  tests 

The  unsteady  tests  covered  a  wide  range  of  Incidences  (-8  to  50  deg) , 
frequencies  (1  to  16  Hz)  and  amplitudes  (1  to  16  deg).  Most  of  the  pres¬ 
sure  and  load  measurements  were  performed  at  80  m/s  wlndspeed,  zero  side¬ 
slip  and  at  harmonic  oscillation,  with  reduced  frequencies  based  on  the 
root  semi-chord  up  to  0.5.  A  survey  of  this  part  of  the  unsteady  testpro- 
gram,  Including  the  runnumbers,  is  shown  In  table  5. 

The  Influence  of  the  wlndspeed  was  studied  by  repeating  a  number  of 
runs  at  30  m/s  (same  wlndspeed  as  In  the  visualization  runs)  and  55  m/s 
(see  tables  6  and  7).  The  frequencies  were  adjusted  at  Che  same  time  to 
correspond  with  the  reduced  frequencies  as  applied  In  the  80  m/s  runs.  At 
80  m/s  also  runs  with  sideslip  were  performed  (see  tables  8  and  9). 

In  all  runs  with  harmonic  oscillation,  the  zeroth  and  first  harmonic 
components  of  the  pressures  and  overall  loads  were  measured,  and  time  re¬ 
cordings  were  made.  A  survey  of  the  runnumbers  for  which  time  histories 
of  the  pressure  distributions  are  presented  Is  given  In  Cable  10.  Table  11 
presents  a  survey  of  the  runnumbers  for  which  time  history  plots  of  the 
overall  loads  at  harmonic  oscillation  are  presented.  In  this  cable  the 
components  of  Che  balance  and  the  corresponding  power  spectral  density 
plots  with  their  plotnumbers  are  Indicated.  In  total,  about  1000  runs 
with  harmonic  oscillation  were  performed. 


In  the  results  of  the  steady  measurements  (see  figure  11)  some  cha¬ 
racteristic  incidence  ranges  can  be  distinguished: 
up  to  9  deg:  attached  (''linear”)  flow 
9  to  19  deg:  fully  developed  vortex  flow 
19  to  36  deg:  vortex  burst  extending  from  trailing  edge 
beyond  36  deg:  vortex  burst  penetrating  the  strake>  almost  fully 
stalled  flow. 

Special  emphasis  was  placed  on  incidences  which  marked  transition 
of  the  flow  characteristics,  or  were  typical  for  the  flow  characteristics 
in  some  incidence  range.  These  actual  Incidences  were  9,  19,  22,  36  and 
42  deg,  which  correspond  to  the  adjusted  incidences  referred  to  in  the 
tables  (aa  ■  8,  18,  22,  38,  46  deg).  At  these  Incidences  a  dense  grid  of 
amplitude  and  frequency  values  was  measured. 

Load  measurements  at  (1-cos)  inputs 

This  part  of  the  test  program  consisted  only  of  unsteady  runs  in 
which  a  (l-cos)  signal  was  used  as  input,  to  simulate  maneuvers.  Time 
history  recordings  were  made  of  the  response  signals  of  the  six  components 
of  the  balance,  using  the  model  motion  signal  as  a  reference. 

The  adjusted  incidence  (oa)  at  the  start  of  the  input  is  varied  from 
8  to  46  deg.  The  duration  of  the  (1-cos)  inputs  is  varied  from  0.083  to 
0.500  s,  the  magnitude  of  the  wing  incidence  variation  was  8,  16,  24  and 
32  deg.  A  survey  of  the  test  program.  Including  the  runnumbers,  is  presen¬ 
ted  in  table  12.  At  the  same  test  conditions,  also  the  first  half  (increas¬ 
ing  incidence)  and  second  half  (decreasing  Incidence)  were  performed  sepa¬ 
rately,  with  some  time  in  between  to  allow  the  flow  to  stabilize. 

Flow  visualization  experiments 

The  intention  of  the  flow  visualization  program  was  to  obtain  infor¬ 
mation  about  the  flow  over  the  upper  side  of  the  model.  This  was  done  at 
the  five  characteristic  incidences,  mentioned  in  section  5.1.2,  both  by 
photographs  and  video  recordings  as  described  in  sections  2.6  and  4,3  at 
a  low  wlndspeed  (30  m/s)  to  obtain  good  smoke  visualization.  A  survey  of 
the  test  conditions  is  given  in  table  13  a  through  c. 

During  the  visualization  tests  the  phase  angle  of  the  pulsation  with 
respect  to  the  model  motion  was  varied  over  a  complete  cycle  in  steps  of 


45  deg.  At  three  chord  positions,  corresponding  with  the  location  of  the 
pressure  transducers,  at  five  incidences,  amplitudes  and  reduced  frequen¬ 
cies  corresponding  with  the  80  m/s  runs,  the  smoke  patterns  In  the  light 
screen  were  registered  on  photofllm  to  determine  the  positions  of  the 
vortex  cores,  as  well  as  on  video,  to  obtain  qualitative  data  of  the  dyna¬ 
mic  flow  phenomenae.  Over  1200  cases  were  completed. 


PRESENTATION  OF  RESULTS 

The  presentation  of  the  test  results  Is  given  In  six  parts  of  the 
report.  The  present  part,  part  1,  contains  some  examples  of  the  test  data 
and  surveys  of  all  testruns.  Some  general  figures,  derived  from  the  test 
are  also  given. 

Parts  II  and  III  present  the  plots  of  the  zeroth  and  first  order  har¬ 
monic  components  of  the  pressure  distributions  for  the  runs  as  presented 
In  the  tables  4  through  9  In  this  part.  In  part  III  the  zeroth  and  first 
order  harmonic  components  of  the  overall  loads,  plotted  vs.  Incidence, 
amplitude  and  frequency  are  also  presented.  The  printed  values  have  also 
been  included  on  microfiche.  All  harmonic  data  are  available  on  magnetic 
computer  tape  In  order  to  enable  easy  data  handling.  The  corresponding 
formats  are  given  In  table  14.  Table  13  gives  an  example  of  a  print  of  a 
steady  test  run  and  table  16  of  an  unsteady  test  run  as  presented  on 
microfiche.  The  test  conditions  have  been  listed  on  the  right-hand  side 
In  the  upper  table.  The  steady  and  unsteady  pressures  of  the  four  sections 
are  given  as  pressure  coefficients  together  with  their  sectional  coordi¬ 
nates.  The  overall  forces  and  moments,  obtained  from  the  balance  are  pre¬ 
sented  as  steady  and  unsteady  aerodynamic  coefficients.  The  sign  defini¬ 
tions  of  the  body-fixed  coordinate  system  are  given  In  figure  1.  The  cor¬ 
responding  pressure  plots  are  shown  In  figures  12  and  13.  The  displace¬ 
ments  derived  from  the  accelerometers  are  presented  both  relative  to 
LVDT  2  and  In  absolute  form.  LVDT  2  is  the  transducer  In  the  model,  where¬ 
as  LVDT  1  Is  the  transducer  in  the  hydraulic  actuator.  When  the  pressure 
transducers  or  the  accelerometers  did  not  operate  properly  the  values  were 
not  presented.  In  the  next  chapter  more  attention  will  be  given  to  these 
harmonic  data. 

Part  IV  presents  the  plots  of  the  time  recordings  of  pressures  and 
overall  loads  for  the  runs  with  harmonic  excitation.  A  survey  of  the  run- 


<1  .*g 


-28- 

numbers  for  which  these  pressure  plots  are  presented  In  part  IV  Is  given 
In  table  10  of  this  part  and  an  example  of  such  a  time  recording  for  the 
four  sections  and  the  eight  phase  angles  Is  presented  In  figure  14. 

Table  11  presents  a  survey  of  the  runnumbers  for  which  the  time  histories 
of  the  overall  loads  are  plotted;  an  example  of  the  overall  loads  vs. 
Incidence,  derived  from  the  time  recordings  of  the  balance  signals  Is 
given  In  figure  15. 

Part  V  presents  a  selection  of  the  power  spectral  density  plots  of 
the  overall  loads  at  harmonic  excitation,  corresponding  with  the  time 
history  plots  of  the  overall  loads  as  presented  In  part  IV  (runnumbers, 
component  and  plotnumbers  are  presented  In  table  11  of  this  part).  The 
responses  of  the  overall  loads  to  (1-cos)  Input  signals  (runnumbers  pre¬ 
sented  In  table  12)  are  also  presented.  An  example  of  such  a  time  history 
Is  given  In  figure  16. 

Part  VI  presents  the  results  of  the  flow  visualization  by  prints  as 
well  as  plots  of  the  variation  of  the  vortex  core  positions  during  a 
cycle  of  the  model  motion  at  various  Incidences,  amplitudes  and  frequen¬ 
cies.  This  part  presents  some  selected  cases  In  tables  17  through  22, 
which  will  be  discussed  further  in  section  7.3.  The  corresponding  photo¬ 
graphs  and  plots  of  the  measured  positions  of  the  vortex  cores  for  the  8 
phase  angles  are  presented  In  figures  31  to  43.  When  measurement  of  the 
vortex  core  positions  was  Impossible  no  values  were  presented  in  prints 
and  plots.  The  photographs  for  the  eight  phase  angles  at  the  aft  light- 
sheet  position,  x/c  >  96.82  Z,  are  also  Included.  Although  it  was  Impos¬ 
sible  to  measure  the  core  positions  In  many  of  these  photos,  they  still 
provide  substantial  qualitative  Information  about  the  flow  phenomena  under 
these  conditions.  Some  plots  (figures  44,  45)  are  also  enclosed  In  this 
part  to  demonstrate  the  Influence  of  amplitude  and  frequency  on  Che  time 
history  of  the  vortex  core  position. 


BRIEF  DISCUSSION  OF  SOME  RESULTS 

In  this  chapter  a  selection  of  wlndtunnel  test  data  Is  discussed.  A 
summary  description  was  also  given  In  reference  12  and  a  preliminary 
analysis  of  harmonic  data  Is  presented  In  reference  13. 
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7 . 1  Overall  loads 

The  variation  of  the  steady  normal  force  coefficient  and  the  moment 
coefficient  with  Incidence  Is  shown  In  figure  11.  The  different  types  of 
flow  are  also  Indicated.  Up  to  about  9  deg,  a  linear  variation  of  CN  and 
Cm  with  Incidence  Is  observed.  Beyond  9  deg  the  slopes  Increase  due  to 
the  development  of  vortex  flows  over  both  the  strake  and  wing.  At  about 
19  deg  Incidence,  vortex  burst  starts  at  the  wing  trailing  edge.  Suction 
at  the  aft  part  of  the  wing  Is  diminished,  while  suction  on  the  front  part 
still  Increases  due  to  the  Increase  of  the  strength  of  the  strake  vortex. 
Beyond  19  deg  Incidence,  the  vortex  burst  moves  upstream,  however,  both 
CN  and  Cm  still  Increase,  though,  at  a  smaller  rate  than  before.  At  about 
36  deg  the  Increase  of  suction  at  the  upstream  part  of  the  wing  due  to 
Increasing  Incidence  and  the  drop  of  suction  at  the  downstream  part  of 
the  wing  due  to  the  growth  of  the  vortex  burst  region,  counterbalance 
eachother  which  leads  to  maximum  values  of  CN  and  Cm.  Beyond  36  deg  the 
flow  over  the  entire  wing  and  strake  collapses  rapidly  Into  an  almost  com¬ 
pletely  separated  flow. 

The  effects  observed  In  the  steady  plots  are  also  reflected  In  the 
unsteady  plots.  This  Is  demonstrated  In  the  figures  17  and  18.  A  very 
distinct  change  In  the  unsteady  pitching  moment  coefficient  Is  observed 
at  about  19  deg  (onset  of  vortex  burst)  and  at  36  deg  where  the  phase 
angle  of  the  unsteady  pitching  moment  changes  about  180  deg.  The  effect 
of  frequency  on  the  first  harmonic  components  of  the  pitching  moment  coef¬ 
ficient  Is  shown  In  figure  19. 

Figure  20  shows  the  real  part  of  the  first  harmonic  of  the  pitching 
moment  coefficient  vs.  frequency  for  several  amplitudes.  At  the  selected 
Incidence  of  a  >  18.9  deg  Re(Cm)  Is  Independent  of  the  frequency  (this 
corresponds  with  the  upper  part  of  figure  19)  but  does  depend  on  the 
amplitude.  This  amplitude  dependence  Is  very  well  organized  at  this  inci¬ 
dence.  Both  frequency  dependence  and  amplitude  dependence  show  up  in  the 
Imaginary  part  of  the  first  harmonic  of  the  pitching  moment  coefficient 
(see  figure  21). 

The  results  of  the  three  force  coefficients  vs.  amplitude  are  shown 
In  the  figures  22  through  24. 

In  figure  22  the  zeroth  harmonic  and  the  real  part  of  the  first  order 
harmonic  component  of  the  normal  force  coefficient  is  presented  versus 
amplitude  at  o  ■  18.9  deg  and  f  3  Hz.  At  this  low  frequency  Re(CN)  de¬ 
creases  with  amplitude  at  the  higher  amplitudes,  as  can  be  expected  from 
the  steady  CN  versus  o  curve,  presented  in  figure  11. 


Figure  23  shows  the  zeroth  harmonic  of  the  side  force  coefficient 
vs.  amplitude  for  3  frequencies  at  a  •  35.9  deg  and  0  >  5.0  deg.  The 
range  used  as  vertical  axis  of  the  graph  Is  less  than  1  Z  of  the  total 
range  of  the  balance  (also  see  section  8.1)»  however,  a  good  distinction 
can  be  made  between  curves  for  various  frequencies.  Figure  24  shows  the 
zeroth  harmonic  of  the  tangential  force  coefficient  vs.  amplitude  for  2 
frequencies  at  18.9  deg  incidence.  The  effect  of  the  frequency  results  In 
a  shift  of  (CT)m  while  the  shape  remains  the  same. 

Figure  25  shows  the  effect  of  sideslip  on  the  unsteady  yawing  moment 
coefficient.  In  the  plot  all  measured  amplitude  frequency  combinations 
for  one  particular  Incidence  are  shown.  The  trends  for  only  one  amplitude 
frequency  combination  are  indicated  In  the  figure.  A  similar  treatment  is 
shown  In  figure  26  for  the  effect  of  sideslip  on  the  first  harmonic  of 
the  rolling  moment  coefficient.  In  the  left  part  of  the  figure  the  fre¬ 
quency  Is  kept  constant  and  the  amplitude  Is  varied,  while  In  the  right 
part  the  amplitude  Is  constant  and  the  frequency  Is  varied. 

Pressure  distribution 

Positioning  of  the  transducers  was  very  appropriate  for  understanding 
the  complex  flow  phenomena  that  were  encountered  In  this  Investigation. 
Thus,  based  on  consideration  of  the  steady  and  unsteady  pressure  distribu¬ 
tion,  It  wa^  possible  to  make  the  flow  field  observations  that  were 
described  In  section  7.1. 

During  the  experiment,  the  zeroth  and  first  harmonic  components  of 
the  pressure  distribution  were  presented.  After  postprocessing,  which  in¬ 
cluded  correction  for  blockage  effects  and  determination  of  the  correct 
incidence,  the  same  presentation  was  used  again.  Figure  12  shows  an 
example  of  the  pressure  coefficients  of  a  steady  test  run  and  figure  13 
shows  the  same  for  an  unsteady  test  run.  In  figure  27  an  overview  is  given 
of  the  development  of  the  steady  pressure  distribution  with  incidence. 

The  chosen  Incidences  are  the  ones  which  are  selected  for  the  flow  visua¬ 
lization  program.  In  section  1  the  presence  of  the  strake  vortex  is  in¬ 
dicated  by  position  of  the  strake  vortex  at  about  45  Z  and  the  wing  vortex 
at  about  80  Z  of  the  local  semi-span.  Up  to  about  19  deg  the  pressures 
grow  with  Incidence.  At  Incidences  greater  than  about  19  degrees  the 
effect  of  the  wing  vortex  on  the  pressure  distribution  decreases.  Beyond 
an  incidence  of  36  deg  the  wing  vortices  have  disappeared  due  to  flow 
separation,  while  beyond  42  deg  the  strake  vortices  have  also  disappeared. 
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At  the  trailing  edge  section  3  shows  the  effect  of  vortex  burst. 
Section  4  shows  the  pressure  distribution  of  the  chordwlse  section  with 
its  characteristic  suction  peak  on  the  nose.  At  incidences  greater  chan 
about  36  deg  trailing  edge  separation  occurs,  while  at  Incidences  greater 
Chan  about  42  deg  flow  separation  on  the  nose  occurs. 

By  combining  plots  of  unsteady  runs,  one  can  easily  get  good  insight 
into  Che  Influence  of  the  different  parameters:  figure  28  shows  the  in¬ 
fluence  of  incidence  and  figure  29  shows  the  influence  of  amplitude  on 
Che  pressure  distribution  at  section  2. 

An  Impression  of  the  contribution  of  the  higher  harmonics  can  be 
obtained  from  the  tape  recordings  by  time  history  plots  (see  figure  14) 
or  spectra.  As  another  example,  figure  30  shows  the  time  history  of  Che 
pressure  distribution  at  section  2  for  the  flow  condition  as  presented  at 
6.98  deg  amplitude  in  figure  29.  The  pressure  signal  at  the  suction  peak, 
which  can  be  associated  with  the  wing  vortex  (80  Z  local  semi-span), 
shows  a  nonlinear  behaviour,  which  is  expressed  in  a  significant  second 
harmonic  component  in  the  spectra. 


Vortex  core  positions 

In  this  section  examples  of  the  time  history  of  the  unsteady  flow 
are  presented  in  figures  31  through  43.  Flow  patterns  have  been  photo¬ 
graphed  for  eight  phase  angles  at  three  spanwlse  sections  and  in  these 
photographs  Che  positions  of  Che  vortex  cores  have  been  determined  as 
shown  in  the  figures. 

At  the  strake  section  (x/cr  40.42  Z)  one  test  condition  was  chosen 
a  -  18.94  deg,  do  *  3.58  deg,  f  ~  1.88  Hz.  The  photographs  of  this  case 
are  shown  in  figure  31  and  the  measured  positions  of  the  cores  are  pre¬ 
sented  in  table  17  and  figure  32.  The  spanwlse  position  of  the  vortex 
core,  depending  on  Che  phase  angle,  is  situated  at  about  70Z,  which  cor¬ 
responds  to  the  steady  pressure  distribution  at  section  1  in  figure  27(3) 
At  these  test  conditions  a  stable  vortex  flow  exists  at  this  section. 

At  the  section  just  behind  Che  kink  (x/cr  ■  65.88  Z)  five  test  condl 
Cions  were  chosen.  Not  only  the  strake  vortex,  but  also  the  wing  vortex 
is  present  in  this  section.  The  situation  in  figures  33  and  34  is  rather 
straight  forward.  At  an  incidence  a  -  18.93  deg  and  an  amplitude  of 
da  «  6.93  deg  we  still  have  stable  vortex  flow  during  the  cycle.  Because 
it  was  attempted  to  optimize  the  smoke  conditions  for  the  wing  vortex, 
the  core  positions  of  the  strake  vortices  were  overexposed  and  could  not 


always  ba  measured  (see  table  18).  The  spanwlse  coordinate  of  the  wing 
vortex  cores  (~  80  Z)  corresponds  with  the  position  of  the  pressure  peak 
In  section  2  In  the  figures  27  and  28.  In  the  figures  35,  36  and  37,  38 
the  Incidence  and  frequency  are  the  same,  but  a  different  amplitude  Is 
used.  Plotting  the  nondlmenslonallzed  z-positlon  of  the  vortex  cores  ver¬ 
sus  phase  angle  (see  figure  44)  shows  that  the  distance  of  the  vortex 
cores  from  the  model  Increases  with  the  amplitude.  The  vortex  flow  at  the 
left  side  breaks  up  at  the  higher  amplitude  at  a  phase  angle  of 
(p  ■  90  deg  (see  figure  37).  In  these  situations  such  smoke  conditions 
were  created  that  the  vortex  core  positions  could  be  determined  for  both 
the  strake  vortex  and  the  wing  vortex  (see  tables  19  and  20).  In  the  test 
conditions  corresponding  to  the  figures  37  through  42  the  same  Incidence 
and  about  the  same  amplitude  were  maintained,  but  3  different  frequencies 
were  used.  In  figure  45  the  nondlmenslonallzed  z-posltlon  is  plotted  ver¬ 
sus  phase  angle.  There  seems  to  be  a  phase  shift  between  the  measurements 
at  1.13  Hz  and  6  Hz,  but  the  measurements  at  3  Hz  mystify  the  picture. 
Also  an  Increasing  frequency  seems  to  postpone  the  vortex  burst  (see  the 
photographs  at  the  phase  angles  la  "  90  and  ip  *•  135  degrees).  The  deter¬ 
mined  positions  are  presented  In  the  tables  20  to  22. 

At  the  trailing  edge  section  (x/cr  ■  96.82  Z)  as  shown  In  figure  43 
the  flow  changes  from  burst  vortices  at  0  degrees  to  fully  separated 
flow  at  a  phase  angle  of  90  degrees,  and  to  stable  vortex  flow  at  270 
degrees  phase  angle. 


FINAL  COMMENTS 

Deformation  and  accuracy  of  the  balance 

For  computer  code  evaluation,  accurate  overall  loads  are  Indispens¬ 
able.  A  balance  Is  used  for  this  purpose,  although  It  will  be  a  weak  ele¬ 
ment  In  the  support  of  the  model.  The  balance  signals  will  also  be  affect¬ 
ed  by  Inertia  loads.  In  this  test  the  resonance  frequencies  of  the  support 
system  (Including  the  balance)  were  far  beyond  the  highest  excitation 
frequency.  Subsequently,  the  Influence  of  these  vibration  modes  on  the 
mudel  motion  were  negligible.  Due  to  the  small  Inertia  properties  of  the 
model  and  the  low  test  frequencies.  Inertia  loads  were  small,  hence  by 
measuring  the  model  motion,  corrections  for  Inertia  loads  were  made  by 
using  wing  model  mass  properties.  In  this  way,  overall  loads  will  be 


obtained  which  are  expected  to  be  nore  accurate  than  airloads  obtained  by 
Integration  of  a  large  number  of  pressures,  measured  at  the  model  surface, 
especially  when  large  pressure  gradients  occur. 

The  maximum  forces  and  moments  which  can  be  measured  with  the  NLR 
635  balance  are  mentioned  In  section  2.3.  Assuming  that  the  balance  can 
measure  accurately  at  less  than  IZ  of  Its  total  range,  the  nondlmenslon- 
allzed  force  and  moment  coefficients  can  be  calculated  corresponding  to 
a  IZ  load.  Comparing  these  values  to  the  measured  values.  It  can  be  con¬ 
cluded  that  the  ranges  of  all  symmetrical  components  conform  to  the  acting 
loads  upon  the  balance.  However  the  loads  of  the  asymmetrical  components 
are  acting  In  the  lower  lOZ  of  the  total  range.  Nevertheless  It  can  be 
said  that  the  accuracy  of  the  measurements  will  be  better  than  IZ  of  their 
total  range. 

Instrumentation 

Just  before  the  test  period,  transducer  3  of  the  42  pressure  trans¬ 
ducers  was  found  unusable  In  both  steady  and  unsteady  measurements.  The 
transducers  6,  27  and  42  were  unusable  In  the  steady  measurements  and  the 
transducers  29  and  37  In  the  unsteady  measurements.  The  measured  values 
of  these  transducers  were  not  presented  In  the  print  out  and  plots.  This 
presentation  was  also  used  with  the  transducers  which  were  unreliable  or 
failed  during  the  tests.  The  transducers  In  pressure  section  1  on  the 
strake  were  highly  susceptible  to  failure  as  a  consequence  of  the  high 
acceleration  levels  encountered  In  this  section. 

Accelerometers  5  and  6  gave  very  poor  signals  before  the  test  period 
and  were  not  presented  in  any  run.  Furthermore  only  the  results  with  acce¬ 
lerometers  1,  7  and  8  seemed  reliable.  The  other  accelerometers  were  not 
reliable  in  all  cases,  which  can  be  seen  In  the  last  column  of  the  LVDT/- 
Acce^erometer  table  In  the  print  out.  In  which  the  calculated  amplitude 
divided  by  the  measured  amplitude  Is  presented. 

Two  of  the  supporting  struts  were  used  as  conduits  for  the  instrumentation 
wiring.  To  minimize  loads  on  the  balance  due  to  the  Instrumentation  wiring 
sufficient  slack  was  provided  between  the  supporting  struts  and  model. 


Backlash  in  vlng  support  bearings 

During  the  wlndtunnel  test  all  signals  could  be  monitored  on  an 
oscilloscope.  By  monitoring  the  LVDT  signal,  an  Increase  In  backlash  in 
the  bearings  of  the  wing  support  could  be  detected.  In  which  case  the 
bearings  could  be  adjusted  and  the  backlash  removed. 


Visualization 

During  the  visualization  at  the  three  spanwlse  sections  a  large 
number  of  photographs  was  taken  with  the  aim  to  measure  the  positions  of 
the  vortex  cores.  This  was  not  possible  In  all  cases;  especially  for  sec¬ 
tion  2  (just  behind  the  kink)  It  was  not  always  possible  to  Inject  the 
right  amount  of  smoke  In  both  vortex  systems,  vnien  more  smoke  was  Injected, 
the  wing  vortices  could  be  photographed  correctly,  but  the  strake  vortices 
were  overexposed;  when  less  smoke  was  Injected  the  wing  vortices  could 
not  be  detected  on  the  photographs.  At  the  third  section,  the  vortices 
were  burst  In  many  of  the  tests  and  only  qualitative  Information  could  be 
obtained  from  Che  photographs. 

The  video  recordings  also  provided  much  qualitative  Information  of 
the  flow,  especially  at  vortex  burst,  because  no  time  averaging  over 
several  light  pulses  was  needed  for  proper  exposure. 


The  effect  of  the  wire  suspension  blocks 

In  one  of  Che  preparatory  tests  as  discussed  In  chapter  3,  measure¬ 
ments  were  carried  out  on  the  model  suspended  by  wires  (see  Ref.  7).  This 
was  done  In  order  to  make  an  estimate  of  the  Interference  of  the  model 
support  system  used  In  the  test  described  In  this  report.  To  be  able  to 
attach  the  wires  In  the  preparatory  test,  special  wire  suspension  blocks 
had  to  be  installed  In  Che  model.  In  figure  46  the  effect  of  Che  wire  sus¬ 
pension  blocks  on  the  lift  and  drag  curves  for  the  model  on  the  strut 
support  Is  presented.  The  graphs  with  and  without  wire  suspension  blocks 
are  presented  next  to  each  other  because  placing  the  two  curves  in  one 
figure  shows  no  difference.  Therefore  the  effect  of  the  wire  suspension 
blocks  can  be  neglected  In  the  zero  lift  correction  as  described  in  sec¬ 
tion  4.2.1. 
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CONCLUDING  REMARKS 


The  aims,  mentioned  in  section  1  have  been  achieved. 

-  A  large  data  base  of  both  overall  forces  and  moments  as  veil  as  pressure 
distributions  has  been  created.  It  Is  expected  that  this  data  base  will 
be  useful  for  the  evaluation  of  advanced  computer  codes  for  the  deter¬ 
mination  of  unsteady  aerodynamic  quantities  for  conditions  at  high 
angles  of  attack. 

-  From  photographs  quantitative  Information  about  the  vortex  core  posi¬ 
tions  Is  made  available  to  better  understand  the  unsteady  vortex  flow. 

-  Video  and  photograph  registrations  made  during  this  experiment  will 
contribute  significantly  to  enlarge  the  physical  Insight  Into  both 
steady  and  unsteady  vortex  flow  phenomena. 
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TABLE  1 

Positions  of  pressure  transducers  and  accelerometers 
(x  100  %) ,  see  figure  5 


PRESSURE  TRANSDUCERS 


no 

x/c 

2y/b 

sec 

no 

x/c 

2y/b 

1 

40.42 

6.81 

III 

29 

96.82 

20.00 

2 

40.42 

20.43 

30 

96.82 

30.00 

3 

40.42 

34.06 

31 

96.82 

40.00 

4 

40.42 

47.68 

32 

96.82 

50.00 

5 

40.42 

54.48 

33 

96.82 

60.00 

6 

40.42 

61.29 

34 

96.82 

70.00 

7 

40.42 

68.10 

35 

96.82 

80.00 

8 

40.42 

74.92 

36 

96.82 

90.00 

9 

40.42 

81.73 

10 

40.42 

88.54 

IV 

37 

5.72 

40.00 

22 

16.61 

40.00 

11 

65.88 

13.11 

38 

29.21 

40.00 

12 

65.88 

26.00 

39 

41.82 

40.00 

13 

65.88 

32.44 

40 

54.42 

40.00 

14 

65.88 

38.89 

41 

67.02 

40.00 

15 

65.88 

42.93 

42 

79.62 

40.00 

16 

65.88 

46.93 

31 

92.22 

40.00 

17 

65.88 

50.98 

18 

65.88 

55.02 

19 

65.88 

59.02 

SECTION 

c 

b/2 

20 

65.88 

63.07 

21 

65.88 

67.07 

I 

785.5 

79.16 

22 

65.88 

71.11 

23 

65.88 

75.56 

II 

785.5 

225.0 

24 

65.88 

80.00 

25 

65.88 

84.44 

III 

785.5 

400.0 

26 

65.88 

88.89 

27 

65.88 

93.33 

IV 

321.38 

400.0 

28 

65.88 

97.78 

ACCELEROMETERS 


no 

x/cr 

2y/b 

1 

82.75 

86.25 

2 

92.94 

86.25 

3 

82.75 

-86.25 

4 

92.94 

-86.25 

5 

92.94 

-37.50 

93-70 

46.72 

21.26 

62.38 


400.0  mm 


2y/b 


.0 
.0 

-.0 

-37.5 


'wlH, 


TABLE  2 

Vibration  modes  with  corresponding  frequencies 
of  Che  Installed  model 


mode 

number 

freq. 

type 

1 

38.66 

roll 

2 

45.36 

roll  +  pitch 

3 

53.03 

pitch 

4 

111.87 

model  bending 

5 

31.97 

yaw 

6 

80.03 

Y-dlsplacemenC  on 
support 

TABLE  4 

Steady  test  program  (B  "  0  deg,  V  •  80  m/s) 


r 


|*ft.t*fa#  't^' 


:; 


TABLE  5  (cont'd) 


s: 


-46- 

TABLE  5  (cont'd) 


V  “  80  m/s 
6-0  dee 


4 

6 

8 

10 

591 

592 

593 

594 

657 

658 

659 

660 

595 

596 

597 

662 

663 

664 

1050 

1051 

599 

666 

•V. 

W 

*Iv 

*1'' 

ri;' 


620 


368 

369 

370 

375 

376 

377 

603 

604 

605 

371 

372 

373 

378 

379 

380 

613 


382 

383 

384 

385 

386 

388 

389 

390 

391 

392 

610 

611 

612 

394 

395 

396 

397 

399 

400 

401 

402 

617 

618 

619 

I 

1 


•  i 
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TABLE  5  (concluded) 


V  ■  80  m/s 
6*0  deg 


a 

a 

nr— 

do 

a 

2  4  6  8  10  12  14  16  18 

54 

3 

423  (T) 

5 

424  (!•) 

8 

641  (!•) 

12 

642  (!•) 

54 

16 

643  (!•) 

-49- 


TABLE  6 

Unsteady  test  program  (6  ~  0  deg,  V 


«  55  m/s) 


V  ■  55  m/s 
g  ■■  0  deg 


I 


“a  1 

f 

8 

2.06 

3.44 

5.50 

8 

11.0 

16 

2.06 

3.44 

5.50 

16 

11.0 

18 

2.06 

3.44 

5.50 

18 

11.0 

20 

2.06 

3.44 

5.50 

20 

11.0 

22 

2.06 

3.44 

5.50 

22 

11,0 

da 

a 

2 

4  £ 

8  10 

12 

14 

16  18 

790 

791 

792 

793 

794 

795 

796 

797 

798 

832 

833 

834 

835 

836 

837 

838 

839 

840 

841 

842 

843 

844 

845 

84o 

849 

859 

847 

848 

850 

851 

852 

853 

854 

855 

858 

856 

857 

860 

861 

862 

863 

864 

865 

866 

867 

868 

883 

TABLE  7  (concluded) 


TABLE  10 

Runnumbcrs  for  which  plots  are  presented  of  tine  histories 
of  the  pressure  distributions 


V  -  80  m/s 
6*0  dee 


l 

0  1 

2 

14  16 

$ 

I 

f 

I 


Pi 


•4*!wwvS*IwwwW'«*?«*»*ni 


fcfft  t*«  i 

4.4  1.'  ir-  $ 

'a* 

-57- 

TABLE 

10  (concluded) 

V  - 

80  m/s 

s  - 

0  deg 

a 

f 

da 

a 

a 

2  3  S’ 

8  10 

12  14 

16 

18 

38 

3 

1044 

1045 

1046 

5 

1047 

1048 

1049 

8 

1050 

1051 

38 

16 

1052 

46 

3 

1053 

1054 

5 

1055 

1056 

8 

1057 

1058 

46 

16 

1059 

t*. 


W.' 


-58- 


TABLE  11a 

Runnuabars  and  corresponding  plotnuabers  of  tine  histories 
and  power  spectral  densities  of  overall  loads  at 
harmonic  oscillation  (V  «•  80  m/s.  6  >  0  deg) 


V  -  80  m/s 
6-0  deg 


3 

N  1108 

115 

n  1109 

17 

- 

129  - 

f 

N  3244 

H 

m  3243 

,*•  **t"i**t<*>*aH^*** 


a 

a 

f 

8 

12 

16 

18 

22 

3 

184 

m 

1125 

186 

m 

1127 

188 

a 

1129 

189 

m 

1130 

5 

202 

m 

1143 

204 

m 

1145 

206 

a 

1147 

207 

D 

1148 

8 

532 

m 

3270 

24 

3 

211 

m 

1152 

213 

m 

1154 

215 

a 

1156 

5 

219 

m 

1160 

221 

m 

1162 

223 

a 

1164 

8 

540 

tn 

4104 

26 

3 

227 

m 

1168 

229 

m 

1170 

5 

235 

- 

237 

a 

1202 

8 

546 

M 

n 

4111 

4110 

28 

3 

242 

m 

1207 

244 

m 

1209 

5 

248 

m 

1213 

250 

a 

1215 

8 

552 

m 

4117 

30 

3 

254 

ID 

1218 

5 

259 

m 

1223 

8 

558 

in 

4123 

32 

3 

262 

m 

1226 

274 

N 

a 

1243 

1242 

5 

264 

m 

1228 

280 

N 

a 

1255 

1254 

8 

564 

m 

4129 

U'« 

5'« 

if* 

Pk 

I 


f 

34 

1 

3 

5 

8 

36 

3 

5 

8 

38 

3 

5 

8 

1 

40 

1 

3 

5 

8 

42 

3 

1 

5 

,  1 

1 

1 

8 

1 

44 

3 

1 

5 

8 

N  2103 
m  2102 
N  4150 
m  4149 


N  2117 
m  2116 
N  2133 
m  2132 
N  4162 
m  4161 


N  2165 
m  2164 
N  2222 
o  2221 
N  5116 
m  5115 


370  - 
377  - 
605 


384 

1 

390 
612  - 


N  5132 


N  2268 
m  2267 
N  3106 
m  3105 


N  3118 
m  3117 
N  3128 
m  3127 

N  5156 
1  5157 
m  5155 


12 


286  ’ 

1 

294  - 


N  1267 
n  1266 


N  2121 

D  2120 

N  2137 
m  2136 


N  2169 
m  2168 
N  2226 
m  2225 


372  - 
379  - 


N  2272 
n  2271 
N  3110 
a  3109 


N  2173 
a  2172 
N  2232 
a  2231 


8 


k  Itk.'ftVk'ltt 


TABLE  11b 

Runnumbers  and  corresponding  plotnumbers  of  time  histories 
and  power  spectral  densities  of  overall  loads  at 
harmonic  oscillation  (V  80  m/s,  8*5  deg) 


V  “  80  m/s 
8*5  deg 


f 

da 

a 

8 

12 

16 

18 

N 

7105 

N 

7108 

8 

3 

800 

1 

7106 

801 

1 

7109 

m 

7104 

m 

7107 

N 

7114 

N 

7116 

5 

803 

1 

7115 

804 

1 

7117 

m 

7113 

N 

7125 

8 

807 

1 

7126 

m 

7124 

N 

8103 

18 

3 

942 

1 

8104 

m 

8102 

N 

8112 

5 

945 

1 

8113 

m 

8111 

N 

8124 

8 

949 

1 

8125 

m 

8123 

N 

7163 

N 

7166 

22 

3 

912 

1 

7164 

913 

1 

7167 

m 

7162 

m 

7165 

N 

7172 

5 

915 

1 

7173 

916 

1 

7201 

m 

7171 

m 

7174 

N 

7209 

8 

919 

1 

7210 

m 

7208 

N 

5241 

N 

5244 

38 

3 

755 

1 

5242 

756 

1 

5245 

ra 

5240 

m 

5243 

N 

5250 

N 

5253 

5 

758 

1 

5251 

759 

1 

5254 

m 

5249 

m 

5252 

N 

5262 

8 

762 

1 

5263 

m 

5261 

N 

6252 

46 

3 

783 

1 

6253 

m 

6251 

N 

6258 

5 

785 

1 

6259 

m 

6257 

N 

6267 

8 

1  788 

1 

6268 

L_ 

m 

6266 

I 


I 

h 
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TABLE  11c 

Runnumbers  and  corresponding  plocnumbers  of  time  histories 
and  power  spectral  densities  of  overall  loads  at 
harmonic  oscillation  (V  ~  80  m/s,  B  -  -5  deg) 


■  80  m/s 
-  -5  deg 


“a 

f 

8 

3 

5 

■ 

8 

18 

3 

5 

8 

22 

3 

5 

8 

38 

3 

5 

8 

46 

3 

5 

8 

8 


N 

810  1 
m 


N 

813  1 

m 


N 

817  1 

m 


N 

932  1 

m 


N 

935  1 

m 


N 

939  1 
m 


N 

922  1 

m 


N 

925  1 

m 


N 

929  1 

m 


N 

765  1 

m 


N 

768  1 


N  7137 
1  7138 
m  7136 


N  7146 
1  7147 
m  7145 


N  7221 
923  1  7222 

m  7220 


N  7230 
926  1  7231 

m  7229 


N  5274 
1  6201 
m  5273 


N  6209 
1  6210 
m  6208 


N 

775  1 

m 


N 

777  1 

m 


N 

780  1 


«*# 


t  •fci'f  .i  •A 


TABLE  12 

Unsteady  test  program  ((1-cos)  Inputs) 
(6-0  deg,  V  -  80  m/s) 


V  -  80  m/i 
6-0  dee 


da 

a 

4  8 

12  16  20 

24 

28  32  36 

3047 

3053 

3058 

3062 

3048 

3054 

3059 

3063 

3049 

3055 

3060 

3064 

3050 

3056 

3061 

3065 

3051 

3057 

3052 

3066 

3067 

3068 

3069 

3013 

3019 

3024 

3028 

3014 

3020 

3025 

3029 

3015 

3021 

3026 

3030 

3016 

3022 

3027 

3031 

3017 

3023 

3018 

3070 

3071 

3072 

3073 

3074 

TABLE  12  (concluded) 


V  ■  80  m/s 
6  ■  0  deg 


32  0.083 


38  0.500 


0.250 
0.200 
38  0.125 


46  0.500 

0.333 
0.250 
0.200 
0.125 
46  0.083 


I 


I 


TABLE  13a 

Visualization  test  program  (photographs)  at  x/cr  »  40.42 


a 

da 

f 

table 

part  VI 

9.98 

4.04 

1.13 

1 

9.87 

8.11 

1.13 

2 

9.28 

16.59 

1.13 

3 

10.01 

3.78 

1.88 

4 

9.91 

7.60 

1.88 

- 

9.38 

15.51 

1.88 

5 

10.00 

3.74 

3.0 

6 

9.88 

7.47 

3.0 

7 

9.42 

15.23 

3.0 

8 

10.00 

3.68 

6.0 

9 

9.88 

7.36 

6.0 

- 

18.96 

3.82 

1.13 

18.92 

7.65 

1.13 

10 

18.78 

13.50 

1.13 

11 

18.94 

3.58 

1.88 

12 

18.93 

7.15 

1.88 

13 

18.79 

12.63 

1.88 

14 

18.97 

3.54 

3.0 

- 

18.92 

7.07 

3.0 

- 

18.83 

12.42 

3.0 

- 

18.98 

3.46 

6,0 

- 

18.93 

6.93 

6.0 

- 

22.45 

3.79 

1,13 

22.41 

7.57 

1.13 

- 

22.29 

15.19 

1.13 

- 

22.46 

3.54 

1.88 

- 

22.42 

7.09 

1.88 

- 

22.29 

14.24 

1.88 

- 

22.44 

3.51 

3.0 

15 

22.42 

6.98 

3.0 

16 

22.28 

14.00 

3.0 

17 

22.50 

3.44 

6.0 

18 

22.42 

6.88 

6.0 

19 

40.42  Z 
30  m 
0  de 


art  VI 


14 

3.73 

1.13 

20 

9 

7.48 

1.13 

21 

13 

15.23 

1,13 

22 

;4 

3.51 

1.88 

- 

7 

7.02 

1.88 

23 

14.26 

1.88 

24 

16 

3.44 

3.0 

25 

7 

6.93 

3.0 

26 

I  *'‘i  *»•  •**  '»1 


-66- 

h 

< 

5 

1 

TABLE  13b 

Visualization  test  program  (photographs)  at  x/cr  - 

■  65.88  J, 

\ 

x/c  -  65.88  Z  1 

V  ■  30  m/s  j 

6-0  deg  S 

a 

da 

f 

table 
part  VI 

9.98 

4.04 

1.13 

28 

9.87 

8.11 

1.13 

29 

9.28 

16.59 

1.13 

30 

10.01 

3.78 

1.88 

31 

9.91 

7.60 

1.88 

32 

9.38 

15.51 

1.88 

33 

10.00 

3.74 

3.0 

34 

9.88 

7.47 

3.0 

35 

9.42 

15.23 

3.0 

36 

10.00 

3.68 

6.0 

37 

9.88 

7.36 

6.0 

38 

18.96 

3.82 

1.13 

39 

18.92 

7.65 

1.13 

40 

18.78 

13.50 

1.13 

41 

18.94 

3.58 

1.88 

42 

18.93 

7.15 

1.88 

43 

18.79 

12.63 

1.88 

44 

18.97 

3.54 

3.0 

45 

18.92 

7.07 

3.0 

46 

18.83 

12.42 

3.0 

47 

18.98 

3.46 

6.0 

48 

18.93 

6.93 

6.0 

49 

22.45 

3.79 

1.13 

50 

22.41 

7.57 

1.13 

51 

22.29 

15.19 

1.13 

52 

22.46 

3.54 

1.88 

53 

22.42 

7.09 

1.88 

54 

22.29 

14.24 

1.88 

55 

22.44 

3.51 

3.0 

56 

22.42 

6.98 

3.0 

57 

22.28 

14.00 

3.0 

58 

22.50 

3.44 

6.0 

59 

22.42 

6.88 

6.0 

60 

35.84 
35.87 
36.01 

35.86 

35.87 
36.02 
35.86 

35.85 

42.32 

42.39 

42.29 

42.39 

42.30 
42.38 

42.31 

42.40 


3.51 

7.02 

14.26 

3.44 

6.93 

14.03 

3.37 

6.79 


m 


I 

li 


r». 

h 
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TABLE  14 

File  organization  on  DELTA  model  tape 


DESCRIPTION 

FORMAT 

RUN. HARM.  ALPHA.  Re(DALPHA).  IM(DALPHA) . 

FREQ.  MACH 

2I5.5F10.5 

VELOC ITY , REDFR . Q , pa . T . BETA , S 

2F10.5.F10.2.4F10.5 

NO.xref .x/xref .yref .y/yref . (Cp)m, 

Re(Cp),  Im(Cp) 

44*(I2.7F10.5./) 

(CN)m,Re(CN) ,Im(CN) . (Cn)m.Re(Cn) .Im(Cn) 

6F(10.5) 

(CY)m.Re(CY) ,Im(CY) , (Cm)m.Re(Cm) .Im(Cm) 

6F(10.5) 

(CT)m,Re(CT) .Im(CT) , (Cl)m.Re(Cl) ,Im(Cl) 

6F(10.5) 

NO.xref .x/xref .yref .y/yref .Re (d) .Im(d) 

9*(I2,6F10.5./) 

N.B.  Improper  values  represented  as:  9999.99 


I 

I 


*  » 
ft 

5*: 

ft 

ft 

a 

rJ 
[< 


1 


Jl>  1>I.  I.  )  A  MUM  :  1  ti> 


..  n  II  It 

I 

0^  £ 

3  a>-*  ft  X 


0^  ^  gi 

i!  V  rt,  N 

■C  C  TJ  X 

'  '“J  K)  O  O  -3 
t  •»  ('J  —  o  »  o 

>  «*  o  ^  o 

^  a 

<  M  tl  I'  II  >i  II 


X  X  O  «  3u. 

<:-  .u  I-  uj  Q 

-  £  ^  JCi  oe 


TABLE  15 

Example  of  a  print  of  a  steady  testrun 


u  tn  -j 

«J  <  Q. 
<  U 
U  £  ^ 

*  »  * 


«  u  » 

X  i  a» 


««««*«*««« 


«««««*«**«« 


«(  «  «  « 


tUu~Uw;jU^L: 


•-  ^  a  •  I 

X  a;  •-< 

ui  O  O  •*  ' 

X  X  X  ' 

X  -1  u 


^  a  ^ 

^  a  a  ^  “j  ^ 

pj  uT  ^  pj  — 


jj  -jT  f*)  —  a  —  *0  -0 

.  «  o  >0  :  s  -0  'O  c  ^ 

*0 


0  ^  c  c 


X  —  ;;  ®  n  “  ' 

x 


-ij  'X  X  ^  o  PJ 

CD'^O  ■C«yrj*«>ts,.*5 


'0  ^  rv  {jj  a 


—  »  ^ 

—  s  ^  -a  t>.  >  > 


'  P'  r*.  —  ~  ^  > 

•3  —  o  —  s  ^  a 


1  -  :i^  -  2  - 


••j  'J  X  -O  -3  - .  V  ^  —  ■^  s  ^  “5  V5  : » 

'  ^  ^4  “?  •■o  ^  '•-•  j-^.  -■;  c  z  -”•  P'  X  72  r  *> 


.LCIitiN  4  i^y/b  O.iOUO  t  =  321.38 


> '  'ftt*  *4  fern*  fc4,* 


I  «*#  »-•  *>t 


TABLE  16 

Example  of  a  print  of  unsteady  testruns 


Q.  a. 


ri  X 
^  ^  Z 

■■3  a—  X 


\  o  c  -o ; 


■J  \ 

•y  ^  -3  r. 


^  'i  >  \ri  « 


^  i-  J  J  -  Ud  ;a 

*  "I  X  -  a. 


j  ®  rj  ^ 
I  —  >  ® 


1  o  33  ^ 
I  a  o  Cd 


i  C  t 
1  U  u 


;r  o  ^  T  ®  4»*) 
o  n N 


«««»«« 


rj  —  rj  ^  ^  if'J  •<0  D 


->  :>  J  W  ;  j  ^  --J  ^  *J  J 


•a 

.c 

•I  -  -  -  - 
-j  a  • 
w 

i  551 


^  ^  D  o.  o  _-5  nj 
>  —  c  j 

■0  •••  <3  ^  ^  ^ 


>  ^  0*  ®  e  ' 

O  >  rj  43  5  .•'J  I 
N  ^  .C  iiT  f*j 


tS  w*)  oT  •-  ® 

3  *«  js  CO  > 

CO  f**  Oj  ''J  *" 


■\  —  ^  i%  “4 

-0  rj  ®  ,0 


i  UT  "O  ^  ^  \  V  S'-' 

-  t  U-)  -C  t .  / 


“  :-v  »'•  3  >  D  \ 

'J  '■'J  “*1  ■*'•  V 


o  ^  ^ 

^  ^  '*  -O 

yl  ^  ro  ^ 


>  iiT  — 


o  >0  -J  ti  —  -J  >  5  -2  I*' ^  ^ 

>  CS  ^  >  *j  '4  —  ^  T  >  v  K  <  '  V  - 

pj  rj  ^  o  -3  o  -c  «*)  T  ^  -0  .1  rj  — j  = 


'  4:  ‘T  o  CK  .’J  -T  -■>  .2  'J  ^  ■'•  “  -'5  -* 

r>^  **5  O  -V  T  —  t’4  >  ■''  o  .:  ‘T  —  '*■  V  "’.’ 

•  -3  'J  C  -J  0  r.  o  >  ^  .2  —  3 


I  ^  ■>  rj  'J  r*3  •■>  J**  oT  r*'  *''  C  ' 


>  *.  gT  *- 

3  O'  “  :5  >  .  — 

^  >  fr  — 


—  0  -3  “4 
_2  ^ 
•3  •=  -  - 


;•  5r  o  Q  X  ^  o  =>  <••  -'  :•  '  "*» 

—  a  7..  .j^  -.  :>  rr  •a-  —  X  •>  —  *2  — 

•r  L"«  w  -.  »  =.  3  r*v 


'4  %  "  -3  3 
0  T  '■■  .:  ' 


>  -.  *4 

•*4  ^  "-  f>  J* 


—  .  O  ^  l"  V 


-3  ^  *2  —  j  ^  r>  o  3  ^  /-.**: 


o3o  ~  = 


-2  ^  „  o  3  -T  .2*  •■  •.  -  O  3  r  "2  o 

' ;  K)  '’•  T  '"  ■)  ■  oT  *»  Jj  —  > 


—  •'J  "  >a*s 


-  2  •v,  r  —  'j  o  C  r** 

^  -  -  oi  <g  <M 


rnw; 


I.*:* 


4 


TABLE 

1  x/c  ■ 

'  b/2  ■ 

1 

40.42 
79.  16 

X  ALPHA  - 

nn  OALPHA  • 

12 

1 

STRAKE 

VORTEX 

1  LEFT  1 

RIGHT 

PHI 

1  2y/b 

22/b  ' 

2y/b 

2r/b 

0 

1 

'  -0.7n 

1 

0.411  1 

0,721 

0.416 

45 

1  -0.710 

0.450  < 

0.71  I 

0.452 

90 

1  -0,719 

0.454  1 

0.691 

0.454 

135 

'  -0.707 

0.443  1 

0.705 

0.446 

180 

'  -0.709 

0.428  1 

0.721 

0.437 

225 

1  -0.716 

0.402  > 

0.723 

0.404 

270 

1  -0.738 

0.381  1 

0.721 

0.386 

315 

!  -0.710 

0.405  1 

0.719 

0.405 

.88  Hz 


UIIN6  yORTEX 


2y/b  2z/b 


2y/b  2z/b 


TABLE  18 

Vortex  core  positions  at  section  2, 
a  ■  18.93  deg,  da  •  6,93  deg,  f  •  6  Hz 


*/c  -  SS,88  X  ALPHA 

b/2  -  225.00  mn  OALPHA 


STRAKE  VORTEX 


18.93  deg  PREO  -  6.00  H: 
5 . 93  deg 


WIN6  VORTEX 


RIGHT 


RI6HT 


K  ' 

0 

! 

1 

-0.782 

0.082 

f 

0.794 

0.092 

P  ' 

45 

1 

1 

-0.782 

0. 1  15 

t 

0.793 

0.127 

i  ' 

90 

!  -0.438 

0.250  1 

0.437 

0.259 

-0.775 

0.153 

( 

0.805 

0.154 

r 

135 

1  -0.442 

0.251  > 

0.442 

0.258 

-0.792 

0.144 

1 

0.814 

0,152 

( 

180 

!  -0.442 

0.229  ! 

0.442 

0.238 

t 

i 

225 

!  -0.444 

0.206  1 

0.443 

0.212 

-0.813 

0.077 

t 

0.820 

0.085 

270 

1  -0.453 

0.180  I 

-0.812 

0.061 

t 

0.824 

0.071 

315 

1 

1 

-0.790 

0.062 

t 

0.331 

0.067 

& 
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TABLE  19 

Vortex  core  positions  at  section  2, 
a  ■  22.45  deg,  da  ■  3.79  deg,  f  *  1.13  Hz 


x/c  -  S5.88  X  ALPHA 
b/2  -  225.00  nn  OALPHA 


STRAKE  VORTEX 


22. 4S  dec  FREQ  *1.13  Hi 
3.79  deg 

- 1  - - - - - - 

!  UIN6  VORTEX 


RIGHT 


-  65.88  X  ALPHA  -  22.41  deg  FREQ  -1.13  Hz  > 

-  225.00  an  OALPHA  -  7.57  deg  i 


STRAKE  VORTEX 


TABLE  21 

Vortex  core  positions  at  section  2, 
a  •  22.42  deg,  da  >  6.98  deg,  f  »  3  Hz 

x/c  •  S5.88  X  ALPHA  -  22.42  deg  FREQ 

b/2  •  225.00  not  OALPHA  ■  6.98  deg 


-3.00  H: 


STRAKE  VORTEX 


WING  VORTEX 


TABLE  22 

Vortex  core  positions  at  section  2  at  22.42  deg 
o  ■  22.42  deg,  da  ■  6.88  deg,  f  •  6  Hz 


1 

TABLE  1 

x/c  -  65.88 

b/2  -  225.00 

X  ALPHA  - 

nw  OALPHA  - 

22.42  deg 
6.88  deg 

FREQ  -  6.00  Hz 

60  1 

STRAKE 

VORTEX 

( 

UIN6 

VORTEX 

1 

LEFT  1 

RIGHT 

1  LEFT 

1  RIGHT 

PHI  I 

2y/b  2z/b  ' 

2y/b  2z/b 

1  2y/b 

2i/b 

'  2y/b  2:/b 

0 

'  -0.469 

0.195 

1 

0.469 

0.205  1 

-0.798 

0.128 

\ 

0.790 

0.133 

45 

'  -0.457 

0.227 

1 

0.457 

0.235  1 

-0.797 

0.150 

\ 

90 

1  -0.449 

0.251 

\ 

0.449 

0.259  1 

I 

135 

1  -0.449 

0.269 

» 

0.449 

0.269  < 

1 

180 

1  -0.448 

0.239 

1 

0.448 

0.250  1 

1 

225 

1  -0.446 

0.210 

1 

0.446 

0.221  1 

( 

270 

'  -0.456 

0.189 

1 

0.456 

0.197  1 

-0.824 

0.085 

1 

315 

'  -0.466 

0.179 

( 

0.466 

0.188  • 

-0.812 

0.  104 

1 

0.81  1 

0.113 

I 


Fig.  6c  Smoke  tube  attachment. 
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m  +dCl: 


MOTION 


am-da 


MODULATION 


N  PULSES 


SHUTTER  OPEN 


SHUTTER  CLOSE 


EXPOSURES  Z  T 


LIGHT  SHEET  ;  PULSE  DURATION  T  INVERSELY  PROPORTIONAL  TO  FREQUENCY 

(  1-8  DEG)  ADJUSTABLE 
PHASE  WITH  RESPECT  TO  MODEL  MOTION 
(0  -  360  DEG)  ADJUSTABLE 


Fig.  10  Principle  of  unsteady  flow  visualization 


LSI* 


Fig.  12  Example  of  a  plot  of  the  pressure  coefficients  of  a  steady  testrun 
(see  also  table  IS) 


Fig.  13  Example  of  a  plot  of  the  pressure  coefficients  of  an  unsteady  testrun 
(see  also  table  16) 
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Exaaple  of  overall  aerodynaaic  loads  va.  incidence,  derived  fron  the 
recordings  at  haraonic  oscillating  aodel 
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Fig.  18  Zeroth  and  first  order  harmonic  components  of  unsteady  pitching  moment 
coefficient  (do  "  3.5  deg,  f  »  3  Ha,  8*0  deg,  V  «  80  m/s) 


Imaginary  part  of  first  order  harmonic  components 
moment  coefficient  vs.  frequency  (o  "  18.9  deg,  8 


Zeroth  order  harmonic  component  of  unsteady  side  force  coefficient  vs 
amplitude  (a  ^  35.9  deg,  8=5  deg,  V  •  80  m/s) 


Fig.  25  The  effect  of  sideslip  on  the  zeroth  end  first  order  harmonic  components 
of  the  unsteady  yawing  moment  coefficient  at  various  frequency  amplitude 
combinations  (a  •>  18.9  deg,  V  ">  80  m/s) 


Fig.  30  Time  history  of  the  pressure  distribution  of  section  2 


a  = 

18.94  deg 

da  = 

3.58  deg 

f  = 

1.88  hz 

Vc  = 

0.4042 

»  A 
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Flf.  32  Tla«  history  of  tho  vortoz  coro  positions  st  section  1 
(o  -  18.04  dS8.  da  -  3.58  dog,  f  >  1.88  Hz) 


Flf.  34  TIm  history  of  tbs  rortex  core  positions  st  ssctlon  2 
(a  >  18.93  ds8.  do  >  6.93  dog,  f  -  6  Hs) 


Fig.  35  Photographs  showing  the  time  history  of  the  flow  at  section  2 
(a  =  22.45  deg,  da  =  3.79  deg,  f  =  1.13  Hz) 
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Fif.  3S  T1«0  history  of  tho  rortox  coro  positiono  st  ooetlon  2 
(a  m  22.41  doc,  da  >  7.87  dof.  f  -  1.13  Rs) 


Fig.  39  Photographs  showing  the  tlae  history  of  the  flow  st  section  2 
(a  >  22.42  deg,  da  s  g.gg  deg,  f  -  3  Hz) 


V08TEXP0S1TI0MPL0T  57 

e«  ■ 

l80-<- 

x  c  -  65.88  •■.  fltPHfl  ■ 

22.42 

4*9 

45-* 

225-0 

b/2  >  225.0  ••  DALPHA- 

6.00 

440 

90-  + 

270- X 

rftCQ  - 

3.00 

Hz 

135-^ 

315-4 

LEFT 


STRQKE 

RIGHT 


se  .4t 

-2v/b 


.4«  .44  .42  .40  .40  .42 


.50 

2y.'b 


LEFT 


UINC 

RIGHT 


+ 

. 

o 

m 

>< 

X 

<  ► 


85  .83  .91 

-2y/t 


.79  .77 


.77  .79 


.85 

2V''b 


Flf.  40  TIm  history  of  tho  rortox  coro  poaltioos  st  soetloa  2 
(a  •  22.42  dog,  da  -  6. Of  dog,  f  ■  3  Rs) 
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Flf.  42  TIm  history  of  ths  Tortoz  coro  positions  st  ssetion  7 
(a  ■  22.42  dsc,  4a  >  6.99  dmg,  f  •  6  Hs) 
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APPENDIX  A 


MODEL  GEOMETRY  AND  DIMENSIONS 


span 

root  chord 
tip  chord 


sweep  angle  strake  leading  edge 
sweep  angle  wing  leading  edge 


area  (wing  +  strake) 


airfoil  of  wing 


aspect  ratio 

oean  aerodynamic  chord 


pitching  axis 


chordwlse 

section 

no. 

2y/b 

c(mm) 

1 

+.30000 

354.715 

1 

-.30000 

354.755 

+.63750 

241.350 

-.63750 

241.410 

+.97500 

128.315 

-.97500 

128.340 

HEHI 

.00000 

785.495 

800.0  mm 
785.495  am 
120.0  am 


76  deg 
40  deg 


.264  in* 
NACA  64A005 


2.422  n 
.320  n 


at  73.29  Zcr 


spanwlse 

section 

no. 

2y/b 

c(mm) 

5 

.12731 

25.000 

6 

.31827 

62.500 

7 

.50923 

100.000 

8 

.87269 

400.000 

9 

.93635 

400.000 

^afloaoagMscwMtt^ 


Note  1:  All  x-coordlnates  are  measured  relative  to  the  leading  edge  of 
the  local  chord  (see  also  figure  A-1) 


TABLE  A1 

Thickness  dlstrlbuClon  of  section  I 


•jtCTION  1 

•3tCTIGN  1 

LtCTlDN  1 

JtCTidN  1 

upptw  siDt; 

LOUbK  SIDb 

UPPtk  SIDE 

LOWLh!  8il;L 

2y/b  3  4-. 3000 

2y/Q  =  •*■.30  0  0 

2y/b  =  -.3000 

2  y  /  D  =  -  .  J  (Hi  u 

c  -  354.715 

c  =  354.715  «« 

c  -  354.755  MM 

c  "  354 , 7' 5 5  MM 

X/C  2/C  X/C  2/C  I  X/C  2/C  X/C  2 / C 


O.OOdOO  -.000;24  .UOUOi  -.OUUb^  O.uOUUO  -.01)047  .uuuio  O.uuOuo 

.00004  0.00000  .00004  -.00092  .00001  0.00000  .00011  -.uUO.o 

.00017  .0003?  .00017  -.00134  .0001  1  .0002b  .uuuM  -  .  u  u  0  jV 

.00025  .00062  .00027  -.00168  .00020  .00062  .OOOlo  -.0u0o2 

.00041  .00092  .00070  -.00227  .00028  .00066  .u0u23  -.uOlib 

.00072  .00124  .00073  -.00244  .00034  .00092  .00042  -OOlbd 

.00100  .00163  .00151  -.00300  .00082  .00145  .00062  -.00213 

.00152  .00197  .00316  -.00379  .00093  .00163  .00109  -.00244 

.00190  .00244  .00488  -.00430  .00163  .00216  .o0l63  -.00279 

.00316  .00307  .00/40  -.00438  .00187  .00244  .0032?  -.00355 

.00489  .00375  .01240  -.00591  .00326  .00313  .00500  -.00407 

.00740  .00450  .02491  -.00818  .00500  .00333  .00751  -.00474 

.01240  .00574  .04993  -.01154  .00751  .00458  .01252  -.00533 

.02491  .00808  .07495  -.01400  .01252  .00531  .02502  -.00827 

.04993  ,01136  .09995  -.01593  .02502  .00806  .02734  -.00874 

.07493  .01380  .14999  -.01900  .05003  .01132  .05003  -.01166 

.09995  .01573  .20002  -.02131  .07504  .01377  .0/504  -.U1409 

.14999  .01879  .25006  -.02302  .10005  .01573  .10005  -.01604 

.20002  .02109  .30009  -.02417  .15009  .01880  .15u09  -.ul?i0 

.25006  .02279  . 350 1 3  - . 02485  .2u011  .02111  .20011  -.02142 

.60009  .02391  .40015  -.02498  .25014  ,02279  ,25014  -.02311 

.35013  .02458  .45019  -.0244o  .30016  .02390  .30016  -.02421 

.40015  .02471  ,50022  -.02344  .35020  .02459  .36021  -.02488 

.45019  .02420  . 55026  -. 021 96  .40022  .02475  .40022  -.025u5 

.50022  .02315  .60028  -.02013  .45025  .02426  ,45025  -.02450 

.55026  .02167  .64236  -.01838  .50027  .02320  ,50027  -.02344 

.60028  .01983  .65616  -.01775  .55030  ,02171  .55031  -.02194 

.65031  .01773  .70035  -.01566  .60033  .01987  ,60033  -.02014 

.70035  ,01538  ,75038  -.01308  .65035  .017/6  .£>5035  -.ul804 

.75038  .01286  .80042  -.01052  .70038  .01543  .70038  -.0156/ 

.80042  .01030  .85044  -.00791  .75041  .01291  ..•5041  -.0  1316 

.85044  .00772  .90048  -.00531  .80044  ,0  1035  .80044  -.OlOtiO 

.90048  .00513  .95051  -.00281  .85046  .00775  .a0u57  -.01053 

.95051  .00262  .99999  -.00044  .90049  .0051&  .85046  -.00798 

1.00000  0.00000  .95052  .00269  ,'»'0u49  -.00533 

1.00000  ,00032  .95052  -,002/o 

I.OUUUO  -.uOOlo 


»>Il  *.J- 


■jtCTIuN  2 

UPPEK  SIDE 

2y/D  =  +.6375 

c  241 . 350  MM 

x/c 

z/c 

U  .  00000 

0.00000 

.00012 

.00062 

. 00025 

.  00087 

. 00027 

.00091 

. 00070 

.00168 

.  00108 

, 00199 

.00155 

. 00244 

.00269 

.00288 

. 00363 

, 00340 

. 00445 

.00360 

. 00694 

. 0U43S 

.01195 

.00559 

.  02447 

.00787 

. 04951 

.01108 

.  0  /454 

, 0 1 353 

.  09957 

.01548 

.  14964 

.01852 

. 19971 

. 02080 

. 24976 

. 02246 

,29983 

. 02360 

. 34991 

. 02428 

. 39998 

. 02443 

.45003 

. 02391 

.50010 

.02289 

.55018 

. 02142 

.60023 

.01960 

.  65030 

.01753 

.70037 

.01523 

.  75042 

.01272 

.80050 

.01013 

. 85057 

.00750 

. 90062 

.00493 

.  95U69 

. 00240 

1  .  00000 

-.00017 

'.itCTiuN  2 


LOUbK  SIDE 
2y /b  =  ♦■ .  6375 


c  ~  241 . 350  Mfi 


0.00000 
.  00004 
. 0002S 
. 00041 
.00106 
.  00114 
. 0026V 
. 00278 
. 00354 
. 0044S 
.  006V4 
.01195 
.  02447 
. 04951 
. 07454 
.  10043 
.14964 
.  19971 
.24976 
.29983 
.34991 
.39998 
.45003 
.50010 
.55018 
.60023 
.65030 
.  70037 
.75042 
.80050 
.85055 
.90062 
. 95069 
. 99985 


0.00000 

-.00093 

-.00128 

-.00168 

-.00230 

-.00244 

-.00327 

-.00340 

-.00379 

-.00398 

-.00476 

-.00605 

-.00833 

-.01162 

-.01403 

-.01601 

-.01898 

-.02132 

-.02300 

-.02418 

-.02480 

-.02490 

-.02438 

-.02337 

-.02188 

-.02003 

-.01782 

-.01560 

-.01305 

-.01044 

-.00783 

-.00526 

-.00276 

-.00041 


JECT iuN  2 


UPPEK  SIDE 


2y/b  =  -.6375 
c  “  241.410 


D.  00000 
.  00008 
. 00031 
.  00054 
. 00124 
.  00137 
. 00298 
.  00313 
. 00474 
. 00721 
.01224 
. 02475 
. 04979 
. 07481 
. 09983 
.14989 
.19995 
.24999 
. 30005 
.35011 
.40017 
.45021 
.50027 
.55033 
.60037 
.65043 
.70049 
.75053 
.80059 
.85065 
.90069 
.95075 
1.00000 


. 00062 
.00000 
. 00168 
.00184 
.00244 
. 00244 
.00325 
.00340 
. 0038? 
.00468 
.00590 
.00822 
.01137 
.01379 
. 01574 
. 01883 
. 02106 
.  02274 
. 02388 
,02450 
. 02456 
.  02405 
.  02307 
.02162 
. 0197» 
.  0 1 769 
. 01541 
.01292 
.01034 
.00777 
. 00524 
.00282 
. 00048 


utLliUN  2 


LOWLK  SIDE 


2y/b  "  -.6375 


c.  “  241.410 


.00027 
00037 
.  UU054 
00058 
.  00097 
00137 
.00180 
.00298 
.  u034a 
.  0  04/4 
.00723 
.01222 
.  02475 
.  04979 
.  07461 
.  08744 
.10308 
.  14989 
.  19995 
.  24999 
.  3uu05 
.35011 
.40017 
.  45021 
.50027 
.55033 
. 60u37 
.65043 
, /0U49 
,75053 
. 80059 
.84124 
. 85 065 
.  86784 
, 89443 
.90069 
.92103 
.  94762 


.00000 
. 00062 
.  0  0  1  U.<4 
. 00073 
. 00 1 08 
.  00205 
. 00244 
. 00309 
.00340 
.  UO j85 
. UU468 
.00592 
. 00824 
.01154 
. 01394 
.01497 
.01613 
.01903 
. 02135 
. 02307 
. 02419 
. 02488 
,  024'? t> 
.  02448 
. 02342 
. 02193 
. 02013 
-.  01804 
-.  01531 

-  .  0  1  3 1  7 
-.  01052 
- .  00833 

-  .  0  0783 

-  ,00690 

-  .  00847 
-.00514 

-  .  0040o 

-  .  00265 
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TABLE  A3 

Thickness  distribution  of  section  3 


jtLTIuN  3 

IjPPfcK  SIDE 

Jv/D  =  +.9/50 

c  120.315  «« 

JECTIUN  3 

LOUEK  SIDE 

2y7b  =  +.9750 

c  =  128.315  «fi 

'-ECTiuN  3 

UPPEK  SIDE 

2v/D  =  -.9750 

c  ”  1 28 .340  MM 

5tLriLlN  3 

LDWEd  SIDE 

V D  =  -  .  9 5  u 

(.  ~  1  28  .  O  <1  l)  MM 

X  /  c 

z/c 

x/c 

z/z 

x/c 

z/c 

x/c 

1/  c 

0  .  uouuo 

0  .  OUOOU 

.  00023 

-.00062 

0 .00000 

.  00062 

. 00023 

0  .  OOOoO 

.00012 

0  .  OOOoO 

.  00027 

-.00094 

.00008 

.00168 

. Oou51 

-  .  0  0  0  6.2 

.00019 

.00070 

.  00047 

-.00168 

.00016  0 

.00000 

. 00055 

- .  0  0  u  y  4 

.  00023 

. 00062 

.  00097 

-.  00195 

.00051 

.00187 

. 00074 

-.00166 

,  ooo^y 

.00094 

.00109 

-.00245 

.00101 

.00245 

. 00129 

-  ,  00230 

.  OOObb 

.00168 

.  00210 

-  .  00339 

.00129 

. 00238 

. uOl 36 

-  .  00245 

. 0009? 

.00163 

.00261 

-.00331 

.00292 

.00308 

.00230 

-  .00339 

.00132 

.00245 

.  0Q323 

-  .  0040  1 

. 00308 

.00339 

. 00292 

-  .  Uu 32 7 

. 002bl 

.00296 

. 00436 

-  ,00401 

.00467 

.00362 

. 00370 

-.0040  1 

.  ou2a4 

. 00339 

. 00686 

-  .  00483 

.00717 

. 00436 

. 004bO 

-  .  U  u  4  0  y 

. 00436 

.  0  03^*7 

.01188 

-.00604 

.01219 

.  00549 

.00717 

-  .  00503 

.  UUddb 

. 00417 

.02439 

-.00822 

. 02470 

. 00783 

,01219 

-  .  U  0a23 

. 0 11 UU 

.00526 

. 04941 

-.01153 

. 04971 

.01103 

.  0247  0 

-  .  008'f2 

,  02439 

, 00764 

. 07443 

-.01383 

. 07472 

. 0 1344 

. 04971 

-  .  0  1  la5 

. 04941 

.01091 

.  09940 

-.01570 

.  09970 

.01543 

.  07472 

-.01406 

.  U7443 

. 01327 

.14944 

- . 01878 

. 14972 

.  01847 

. 09970 

-  .  0  1 1>  U  7 

.  09940 

.01531 

.  19947 

-.02112 

. 19974 

.  02077 

.  14972 

-  .  0  1  9  1 

.  14944 

.01839 

.24950 

-.02272 

.24977 

. 02244 

. 1 9974 

-.  02151 

.  1994? 

.  02073 

.29954 

-.02389 

. 29979 

.  02361 

, 24977 

02326 

.  24950 

. 02241 

. 34957 

-.02451 

.34981 

.  02431 

.29979 

-  ,  02443 

. 29954 

,  J2369 

.39960 

-.02459 

.39984 

.  02454 

. 34981 

-  ,  02505 

. o4957 

. 02439 

. 44964 

-. 02404 

.39984 

, 02447 

. j9y84 

-.02517 

. 39960 

. 02455 

. 49967 

-.02295 

. 44986 

.  02392 

. 44986 

-  ,  02470 

. oVYbO 

.  02459 

.54970 

-.02143 

. 49988 

. 0228/ 

.  49938 

-  .  0  2  3  o  1 

. 44956 

.  02412 

. 59973 

-.01952 

.54991 

.  02131 

. 54991 

-  .  02221 

.  49967 

. 02307 

.64973 

- .01738 

.59993 

.01948 

, 29993 

-  .  0204 1 

. 54970 

. 02159 

.69820 

-.01500 

. 64991 

.  01738 

.  64yy 1 

-. 018j1 

.39973 

.  01980 

.69976 

-.01492 

.69994 

.01480 

.  69y  '4 

-.'>1570 

. 64973 

.  0  1  777 

, 74976 

-.01243 

. 74996 

.01223 

,  7499i> 

- ,  01321 

.69976 

. 01531 

. 79983 

-.00986 

. 79998 

. 00978 

. 79998 

-  ,  0  104  4 

. 749b0 

.01270 

.84986 

-  00721 

.85001 

.0  0729 

. 85001 

-  .  0  0  775 

. 799B3 

.01009 

.89989 

- ,00460 

.90003 

. 00483 

. 89999 

-  ,  0  051  4 

. B49b6 

.00740 

,  94993 

- , 00218 

.95005 

.00242 

,  95005 

-  .  002B4 

.  139989 

. 00483 

,  99852 

-.00019 

1.00000 

.  00312 

. 99856 

-  ,  UU09U 

.94993 

.00245 

.  99844 

.  00051 

1.00000 

.  00312 

TABLE  A4 

Thickness  distribution  of  section  4 


SECTION  4 

upper  side 

lower  side 

2y/b  •  0. 

0 

2y/b  -  0.0 

cr  *  785. 

495  ma 

cr  ■  785.495  nm 

x/cr 

r/cr 

x/cr  z/cr 

.57547 

.02798 

.57685  -.02806 

.72125 

.02915 

.73357 

.02798 

.82962 

.02799 

.81258  -.02805 

TABLE  AS 

Thickness  distributions  of  sections  5  through  9 


UPPER  SIDE 

LOWER 

SIDE 

UPPER 

SIDE 

LOWER 

SIDE 

SECTION  5 
x/cr  -  .12731 
b/2  *  25  mm 

2y/b  2z/b 

2y/b 

2z/b 

2y/b 

2z/b 

2y/b 

2z/b 

.93420 

.01140 

.01100 

-.93420 

BSI 

.54220 

.09060 

.54220 

-.08800 

-.54220 

.08980 

-.54220 

-.08840 

.22840 

.15360 

.22840 

-.15140 

-.22840 

.15280 

-.22840 

-.15140 

SECTION  6 
x/cr  -  .31827 
b/2  *  62 •  5  iniD 

2y/b  2z/b 

2y/b 

2z/b 

2y/b 

2z/b 

2y/b 

2z/b 

.97752 

.00320 

.97760 

-.00416 

-.97752 

-.97752 

-.00416 

.57048 

.08504 

.57048 

-.08608 

-.57048 

-.57048 

-.08608 

.25032 

.14944 

.25032 

-.15016 

-.25032 

.14992 

-.25032 

-.15016 

SECTION  7 
x/cr  -  .50923 
b/2  -  100  mm 

2y/b  2z/b 

2y/b 

2z/b 

2y/b 

2z/b 

2y/b 

2r/b 

.98155 

.00610 

.98155 

-.00720 

-.98155 

.00605 

-.98155 

-.00690 

.53045 

.09340 

.53045 

-.09470 

-.53045 

.09360 

-.53045 

-.09435 

.09770 

.18055 

.09770 

-.18180 

-.09770 

.18115 

-.09770 

-.18135 

SECTION  8 
x/cr  -  .87269 
b/2  ■  400  mm 

2y/b  2z/b 

2y/b 

2z/b 

2y/b 

2z/b 

2y/b 

2z/b 

.15000 

.00241 

.15000 

-.00244 

-.15000 

.00241 

-.15000 

-.00242 

.07500 

.00340 

.07500 

-.00343 

-.07500 

.00340 

-.07500 

-.00342 

SECTION  9 
x/cr  -  .93635 
b/2  ■  400  mm 

2y/b  2z/b 

2y/b 

2z/b 

2y/b 

2z/b 

2y/b 

2z/b 

.15000 

.00120 

.15000 

-.00122 

-.15000 

-.00120 

-.15000 

-.00122 

.07500 

.00169 

.07500 

-.00172 

-.07500 

-.00170 

-.07500 

-.00171 

SECTION  5 
SECTION  6 
SECTION  7 
SECTION  8 
SECTION  9 


SECTION  4 
SECTION  3 
SECTION  2 
SECTION  1 


Fig.  A1  Coordinate  reference  system  and  section  positions 
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APPENDIX  B 

NON-AERODYNAMIC  LOADS  ON  THE  MODEL  BALANCE 


CONTENTS 

1  LOADS  ON  THE  MODEL  BALANCE 

2  MASS  AND  INERTIA  PROPERTIES 

3  VARIOUS  REMARKS 

1  Figure 


LIST 

OF  SYMBOLS 

alf 

wing  Incidence 

(deg, rad) 

alf 

wing  angular  velocity 

(rad/s**2) 

•  • 

alf 

wing  angular  acceleration 

(kgm) 

C 

mo^ 

model  Constant  (H*dx  )  (see  section  2) 

c.g. 

(kgm) 

C 

moj 

model  Constant  (J  -M*dx  *dx.  )  (see  section  2) 

y  c.g.  ^.c 

(kgm**^) 

dx 

distance  In  wing  reference  plane  (see  figure  B.l) 

(m) 

f 

frequency 

(Hz) 

g 

gravitational  acceleratlou 

(m/8**2) 

J 

y 

moment  of  Inertia  about  axis  parallel  to  y-axls 

(kgm**2) 

1 

wing  rolling  moment  (see  figure  B.l) 

(Nm) 

m 

wind  pitching  moment  (see  figure  B.l) 
ref.  axis  through  balance  center 

(Nm) 

M 

mass  (see  section  2) 

(kg) 

N 

wing  normal  force  (see  figure  B.l) 

(Nm) 

n 

wing  yawing  moment  (see  figure  B.l) 

(Nm) 

T 

wing  tangential  force  (see  figure  B.l) 

(N) 

t 

time 

(s) 

X 

chordwlse  coordinate  In  wing  reference  plane; 
apex;  x*0  (see  figure  B.l) 

(m) 

Y 

wing  force  in  y-dlrection  (see  figure  B.l) 

(N) 

y 

spanwlse  coordinate  in  wing  reference  plane 
(see  figure  B.l) 

(m) 

z 

coordinate  in  plane  of  synnetry  normal  to  wing 

(m) 
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SUBSCRIPTS 

0  zeroth  order  hamonic  conponents 

1  first  order  hamonic  conponents 

b  balance 

b. c.  balance  center 

c. g.  center  of  gravity 

g  gravity 

1  Inertia 

mo  model 

r.a.  rotation  axis 


1  LOADS  ON  THE  MODEL  BALANCE 

The  balance  In  the  model  measures: 

-  aerodynamic  loads 

-  Inertial  loads 

-  gravitational  loads 

-  wireloads  caused  by  the  electrical  wires  of  the  miniature  pressure 
transducers  and  accelerometers  leading  from  the  model  to  the  support 
system. 

From  wlnd-off  tests  It  was  concluded  that  the  wireloads  were  ex¬ 
tremely  small.  Therefore  no  correction  for  these  loads  was  applied. 

Gravitational  loads  on  the  balance  are: 

N  -  -M  *  g  *  cos(alf)  (1) 

S 

T  -  +M  *  g  *  sln(alf)  (2) 

S 

m  ■>  -fM  *  g  *  (dx,  -dx  )  *  co8(alf)  (3) 

g  o.c.  c.g. 

Inertial  loads  on  the  balance  are: 


•  c  •  • 

N  -  -M  *  dx  *  alf  -  -C  *  alf  (4) 

1  c.g.  moj 
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T  -  -M  *  dx  *  (alf)* 
1  c.g. 


-  -C  *(alf)* 
no^ 


m.  ■  m.  +  N.  *  (dx  -dx.  ) 
^  ^c.g.  ^  ‘=-«- 


I  *  alf  -  M  *  dx^  ^  *  alf  *  (dx  -dx^  ) 

^c.g.  ‘^•8* 


-  {-J 


-  M  *  dx*  }  *  alf  +  M  *  dx  *  dx.  *  alf 

C»g«  C.g.  • C • 


-J  *  alf  +  M  *  dx  *  dx^  *  alf 

y _  c.g.  T>.c. 


-C  *  alf 
m02 


When  the  wing  performs  a  harmonic  pitching  oscillation*  described  by 
alf  •  *  si"**  (2.pi.f.t), 

the  zeroth  and  first  order  harmonic  components  of  the  inertial  loads 
become : 

N,  -  C  *  alf,  *  ((2.pl.f)**2]  (7 

ll  moj  1 


^4  -  -i  *  C  *  [alf,.2.pi.f]**2 

Iq  moj  1 


m,  -  C  *  alf,  *  [(2.pl.f)**2] 

ij  mo2  1 


MASS  AND  INERTIA  PROPERTIES 

This  section  describes  how  the  constants  M,  C  ,  C  ,  dx  , 

moj  mo2  c.g. 

dXj^  which  show  up  In  the  formulae  In  section  1,  are  determined  for 
this  test  setup. 

Before  the  start  of  the  wlndtunnel  test,  a  very  small  mass  was  added 
to  one  of  the  wing  tips,  to  position  the  center  of  gravity  In  the  plane 
of  symmetry  of  the  model. 

When  the  model  Is  mounted  on  the  balance,  a  part  of  the  balance  may  be 
considered  to  move  ss  an  Integral  part  of  the  model. 


If 


ll 

fir 


I 


i:;*: 


iji' 

•K' 

‘I't 

I 

k 

n 

if*! 

Vi' 

•y; 


w 

w 

n 

if. 

*1*' 


I 

I 

[»%'• 


rJi 


That  part  yields  additional  gravitational  and  Inertial  loads  on  the 

balance  itself.  Therefore,  In  the  determination  of  the  gravitational  and 

Inertial  loads,  the  mass  (M)  and  moment  of  Inertia  (C  )  and  the  loca- 

mo2 

tlon  of  the  center  of  gravity  (dx  ) ,  of  the  wing  and  the  part  of  the 

c.g. 

balance  clamped  to  the  model,  must  be  used.  These  values  are  determined 
as  follows: 

-  the  mass  Is  determined  by  positioning  the  model  on  the 
balance  at  different  angles  of  attack  and  applying  the 
formulae  (1)  and  (2)  - >  M 


by  knowing  the  mass  of  the  model  only,  also  the  mass  of 
the  part  of  the  balance  Is  known.  Form  the  construction 
of  that  part  Its  location  of  the  c.g.  can  be  determined. 
In  combination  with  the  c.g.  and  the  mass  of  the  wing 
only,  the  c.g.  corresponding  to  the  mass  M,  Is  determined 


- >  dx 


-  dx  *  M  •  - >  C 

c.g.  n 

-  by  oscillating  the  model  (at  wlnd-off) ,  the  constant 

C  Is  determined  — c 


the  location  of  the  balance  center  Is  determined  by  the 
construction  of  the  balance  and  the  positioning  of  the 
balance  with  respect  to  the  model 


- >  d: 


*b.c. 


VARIOUS  REMARKS 


1.  In  the  design  of  the  model.  It  was  decided  to  place  the  balance  center 
a  little  forward  of  the  rotation  axis,  providing  a  more  smooth  con¬ 
touring  of  the  thicker  part  of  the  model,  containing  the  balance,  on 
to  the  thin  trailing  edge.  Consequently,  also  the  center  of  gravity 
moved  forward  and  so,  both  the  b.c.  and  c.g.  were  not  located  In  the 
rotation  axis  and  the  corrections  (3),  (4)  and  (5)  (see  section  1) 
became  necessary. 

2.  The  right  way  to  determine  the  moment  of  Inertia  C  would  have  been 

moj 

vibrating  the  model  In  vacuum,  measuring  the  output  m^  of  the  balance 
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and  alf  and  application  of  formula  (6).  However,  the  measurement  was 
done  In  the  wlndtunnel  at  wlnd-off  conditions.  The  effect  of  still 
air  reactions  was  estimated  to  be  small. 

3.  Because  the  b.c.  was  not  located  In  the  rotation  axis,  deformation  of 
the  balance  would  cause  a  change  In  the  location  of  effective  rota¬ 
tion  axis.  This  effect  was  assumed  to  be  small. 

4.  The  maximum  absolute  value  of  the  Inertial  and  gravitational  loads 
are : 


N 

g 

-  70 

(alf  -  0  deg) 

T 

g 

-  50 

(alf  -  50  deg) 

m 

g 

-  0.35 

(alf  -  0  deg) 

-  75 

(aif  -  300) 

^^1 

-  75 

(alfj 

-  1.75  deg. 

^1 

-  16 

e 

(alf 

-  8) 

^1 

^0 

-  8 

(alfj 

-  15  deg,  f 

”l 

-  36 

•  • 

(alf 

-  300) 

"1. 

-  36 

(alfj 

-  1.75  deg. 

5.  The  ratio  between  aerodynamic  loads  and  non-aerodynamlc  loads  depends 
on  the  Incidence,  angular  velocity  and  angular  acceleration.  For  the 
major  part  of  the  test  runs,  the  non-aerodynamlc  loads  are  small 
(less  than  10  Z)  as  compared  to  the  aerodynamic  loads. 

6.  Only  the  symmetrical  components  are  corrected  for  Inertia  effects. 
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/  HORIZONTAL  LINE  IN 

. PLANE OF  SYMMETRY 

^  OF  THE  MODEL. 


7 


ROTATION  AXIS 


a.d.di' 

WING  REFERENCE  PLANE 


Bl  Body-fixed  coordinate  eystee 

Note:  *  y-axla  through  balance  center 

*  wing  reference  plane  la  Indicated  In  figure  5,  Part  I 


APPENDIX  C 

UPDATED  VALDES  OF  THE  PITCHING  MOMENT  COEFFICIENT 
OF  THE  MODEL  SUSPENDED  IN  VIRES 

In  report  NLR  TR  86047  C  (reference  7)  concerning  "Force  neasure- 
nents  In  a  low  speed  vlndtunnel  on  a  nodel  of  the  straked  wing;  suspended 
In  wires",  the  values  of  the  pitching  aoaent  coefficient  were  affected 
by  the  drag  of  the  wires,  which  yield  an  additional  pitching  moment. 
Tables  Cl  through  C7  and  figure  Cl  present  the  correct  pitching  moment 
coefficients  derived  from  only  the  force  In  the  wire  attached  to  the 
front  part  of  the  model. 
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